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FOREWORD 


This report presents the results of a theoretical study of hypervelocity 
impact into semi-infinite targets. The work has been sponsored by the Lewis 
Research Center of the National Aeronautics and Space Administration under 
Contract No. NAS 3-2121, in support of their Research Program to develop a 
structurally adequate radiator for use on space vehicles. In order to provide 
the desired support, diverse aspects of the meteoroid impact problem have 
been examined. This report collects all pertinent data obtained under the 
contract into one reference work which will be useful, as such, for NASA's 
program and also for background for further research on hypervelocity impact 
phenomena. 

Technical monitoring was provided by Mr. Robert J. Denington and 
Mr. James J. Kramer, to whom the authors are very grateful for numerous 
helpful suggestions. 

The approach on which this study is based was suggested originally by 
Dr. Franklin K. Moore of this Laboratory, and was later elaborated by Dr. 
Walter E. Gibson. The authors have continued to benefit from frequent dis- 
cussions with Drs. Moore and Gibson, and in addition, valuable contributions 
have been made by Mr . Howard A. Scheetz and Dr. NormanS. Eiss. 

Special thanks are due to Mr. Harold M. Rosenbaum, who very ably 
handled the programming of various calculations for an IBM 704, and to 
Miss Marcia J. Williams, Miss Sarah J. Geraci, Mr. Girard A. Simons, 
and Mr. John R. Moselle, who prepared the figures and carried out the 
associated calculations. 

The analyses presented below embrace both fluid-mechanical and solid- 
mechanical considerations. The former are the work of the first author, while 
the latter are due to the second. 
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ABSTRACT 


(6 Ml 

This report presents an analytical treatment of the process of crater 
formation in semi- infinite targets as a result of impact by hypervelocity 
projectiles. The results are achieved by using blast-wave theory to treat 
the fluid-mechanical phase of the target response, together with a simple 
means of accounting for the target strength in the later stages of the cratering 
process. The crater- size formula finally deduced is quite simple. To enable 
the interested reader to apply it without following a lengthy derivation, the 
formula and its application are presented in a separate section. 

The theory in its present form is found to agree reasonably well with 
relatively low- speed experiments, and the agreement can be expected to 
improve at higher impact speeds. The simple crater-size formula reveals 
that the kinetic energy of the impacting particle is a controlling parameter 
and that the dynamic strength of the target is the factor most effective in 
limiting penetration. 

Careful attention is paid to energy and momentum conservation, and 
the momentum of the incident particle is found to play a relatively minor 
role. The attendant physical and mathematical reasons are discussed, 
including their implications on experiment. 

Requests for copies of this report should be referred to: 

National Aeronautics and Space Administration 
Office of Scientific and Technical Information 
Washington 25, D. C. 

Attn: AFSS-A 
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INTRODUCTION 


The purpose of this research has been to achieve an analytical 
description of the phenomena brought into play when a high-speed projectile 
strikes a semi-infinite target. The extremely high pressures generated 
under the conditions of high-speed impact suggest a fluid-mechanical 
model for describing the motion of the impacted medium. Numerical 
solutions of the resulting equations have been published for specific cases. 
This report describes analytic results obtained by applying blast-wave 
theory to solve these fluid-mechanical equations. The spirit of the approach 
is to simplify the analysis wherever possible by making certain approxima- 
tions to the true physical situation. The aim in doing so is to achieve 
generality and simplicity of the results, at the expense of some exactness 
in describing the details of the solution. 

The report is divided, essentially, into two parts. The first of these. 
Part I, describes the blast-wave solution for the early, high-pressure phase 
of the motion. The second, Part II, discusses the solid-mechanical aspects 
of the problem, particularly from the point of view of the target strength. 

Part III summarizes the theory, and illustrates the use of the crater- size 
formula by working out several examples. To facilitate quick application, 
this section is so written that it may be read independently of the remainder 
of the report. 

In Part I, two types of solution are described. In one, which permits 
spatial variations in two directions, it is possible to conserve both the energy 
and momentum of the system. The other solution, which allows only for 
variations in the radial direction, and conserves only the total energy, is a 
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modification of the classical Taylor solution for an instantaneous point 
release of energy. Approximate solutions of the former type are found to 
be very close to the vastly- simpler Taylor solution in all important respects. 
The analytical results that lead to this conclusion are described, and the 
physical reasons that underlie this phenomenon are discussed. The 
significance of this finding is to indicate why the simpler Taylor solution 
may be used. 

Part II is devoted generally to a study of the limits of the blast-wave 
analysis, and describes some of the detailed phenomena that may be expected 
during the early stages of the impact process. In addition, particular 
attention is paid to the later stages of crater formation, during which the 
plastic and elastic response of the target must be considered. An important 
part of this work is to identify the material- strength level at which the blast- 
wave solution is to be terminated. The meaning of this as a crater -formation 
criterion is discussed. 

The concluding remarks summarize the advances that have been made, 
and indicate the areas where improvements in the theory can be expected. 


RM-1 655-M-4 


2 



I. BLAST -WAVE SOLUTION 


This portion of the report presents the solutions obtained for the target 
response in the fluid regime. Section A, below, describes the fluid-mechanical 
model, listing the several approximations that are fundamental to a blast-wave 
solution. The second and third sections then present two different types of 
solution of the fluid- dynamical equations. 

A. Fluid-Mechanical Model 

To derive an analytic description of hyper velocity impact, two steps are 
taken. The first is to treat the target material as a compressible fluid, while 
the second is to simplify the resulting equations in such a manner that a simple 
solution is possible. Similarly, below, the basic equations appropriate to the 
compressible-fluid approximation are first discussed. Then described are 
the key mathematical simplifications used, and the restrictions which they 
imply. The last two sub-sections summarize the formulation of the problem. 

1. Basic Equations. When a particle strikes a target surface at high 
speed, large amounts of energy and momentum are quickly deposited over a 
very small portion of the surface. This release drives a strong shock wave 
into the target, generating extremely large pressures, typically measured in 
megabars. Because these pressures are so large compared with the mate- 
rial strength, even at high strain rates, one is led to the approximation 
that the impacted medium behaves like an inviscid, compressible fluid. In 
actual fact, the justification for such an approximation is not provided by the 
magnitude of the pressures themselves, but must come from a consideration 
of their gradients. Consider a small mass element 
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Art 





'd O' d 

The net force acting in the x-direction is proportional to r . 

Thus the neglect of resistance to shear deformation requires !£/ '|£»/ ■ 

/ ^ 

To replace this by the simpler statement above is to assume that rates of 
change in the two perpendicular directions are of the same order, and that 
the proper orders of magnitude to use for p and <y are the impact pressure 
and material strength. There appears to be no reason for doubting either of 
these assumptions in the early stages of the impact process. 

Thus the problem of determining the response of the target material 
becomes essentially that of solving the fluid-mechanical equations expressing 
the conservation of mass, momentum, and energy, together with the equation 
of state of the medium 


Conservation of mass 


Dp ^ 

+ f =0 


(i) 


Conservation of momentum 

D ' 




Conservation of Energy 

Dt r th- 


0 or 


Equation of State 


Ds 

Dt 


- 0 


( 2 ) 


(3) 


e = f far) 


(4) 
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in which ~t is the time and V the gradient operator. It should be noted 
that the assumption of an inviscid fluid has been made by setting the right- 
hand side of Eq. (2) equal to zero. If shearing forces were to affect the motion, 
they would have to be added to this equation. Consistent with this approximation, 
energy changes arising from viscous dissipation and heat conduction are omitted 
from the energy equation. In addition, energy changes due to radiation and 
chemical change are also neglected. Thus, the conservation of energy simply 
states that, for each element of mass, changes of internal energy, aCe , are 
balanced by changes in the flow-work term, -ft* ^ • Alternatively, 

this may be expressed by stating that the entropy of a given mass element does 
not change, after it has been processed by the shock. 

Finally, it should be noted that the use of an equation of state implies 
an assumption of thermodynamic equilibrium. 

2. Symmetry Assumptions . All of the analyses of this report assume 
that the impact occurs at normal incidence, and the ensuing motion is taken 
to be symmetric about an axis which is normal to the original target surface, 
and passes through the impact point. For such an axisymmetric flow, the 
scalar forms of the equations of motion, in spherical coordinates, are 
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Continuity 


±£ + a !d t «!£ +/0 (§tt + j r ^ + ^ + ^ cot ©)=o 

3t dr r 3 Q < \dr r d@ r r / 

Radial component of the momentum equation 

l £ i .0 


<9a , , 9a j_ 9m 

9t * dr r dQ 


ttr 


r p dr 
Angular component of the momentum equation 


'd'ur vr_ 

9t 9r r 90 r 




— = 0 

r 90 u 


Conservation of Energy 


de , ,, o>e , mt 9e 
at 9r r dQ 


-*l 2£. + u Ir + ^ 4e) = 


9t 


9r 


90 


0 


Equation of State 


= F 


( 6 ) 


(7) 


(8) 


(9) 


( 10 ) 


Here U. and mt denote the velocity components in the r- and Q- 
directions, respectively. Equations (6) - (10) are five relations for the quan- 
tities -p , jO , U , aaT' , and e . One can also work with the entropy, rather 
than the internal energy, in which case the last two equations are replaced by 


3s + a 9s + — _ n 

et + dr r 30 ' 0 




( 11 ) 

( 12 ) 


3. Boundary Conditions. The boundary conditions that apply at the shock 
wave, i. e. j at r = t) > state that the discontinuities in velocity, pressure, 

density, etc. across the wave are given by the Rankine-Hugoniot relations. 

For a shock advancing into a medium at rest, these are 
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U, 




= o 




/°o > A > ft o 


Conservation of mass 


C 5 . “s = /*, (“ 5 - «,) 

Conservation of momentum 


- P, + /°, 


(13) 


(14) 


or, using (13) 


* i°° «, 


(14a) 


Conservation of Energy 


e 


i 







(15) 


These three equations, together with the equation of State, are sufficient to 
define the four quantities U t , f, , fa, , and e ( as functions of U $ , for 
given initial values fa o , p , and e Q . This locus of states which may 
be reached across a normal shock wave is called the Hugoniot of the medium. 

The analysis of this report assumes that the shock wave is always hemi- 
spherical in shape as it advances into the target. This assumption is based on 
observations of shock shape in Lucite, in metals, ^ and in wax, 

all under hypervelocity impact conditions. Further evidence for the assumption 
comes from fact that the craters formed are very nearly hemispherical in shape 
at high impact speed. 
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At this point, then, the fluid-mechanical problem that is posed is the 
solution of Eqs. (6) - (10), describing the motion of an inviscid, compressible 
fluid behind a hemispherical shock wave that advances into a semi-infinite 
target 


The motion must be such that the boundary conditions (13) - (15) are satisfied 
at the shock, while along the surfaces AB and CD (whose location is unknown) 
the pressure and material density must vanish. 

The solution of such a boundary- value problem is an exceedingly for- 
midable task; the coupled, nonlinear partial differential equations involve 
three independent variables. Approximate numerical solutions of these equations 
have been presented by Bjork, ^ for some specific cases. The objective of the 
present work is to achieve approximate analytic solutions which display, in a 
simple but realistic way, the influences of the various physical parameters. 

4. Similarity Assumption. The most important approximation made, 
in order to achieve the goal of an analytic solution, is that the flow is self- 
similar; i. e. , the distributions of the various physical quantities (such as 
pressure, density, etc. ) at each instant, are taken to be the same when 
viewed on a scale defined by the shock radius at that instant. Thus each 
quantity, instead of depending separately on the time and on the distance r 
from the impact point, is assumed to be a function of the combination 
ir ^ . This reduction of the number of independent variables 
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constitutes a significant simplification in the differential equations that must 
be solved. 


The mathematical expression of this assumption is to introduce the 
similarity variable 

^ " ejit) (i6) 

and to redefine the velocity components, pressure, density, and internal 
energy by the dimensionless functions 

u ( r > o, t) = R s <p(r\, 6) -p[r, Q,f) = p 0 R*-f (rj, q) 

«r(r,e,t)= R s co(ri,e) P(r,e,{) = p 0 (1?) 

e(r,6,t)= 


When these relations are substituted into Eqs. (6) - (9), and derivatives with 
respect to r and t replaced in terms of derivatives with respect to rj , 
one finds that all explicit time dependence disappears from the differential 
equations if one chooses 


R, = At* 

When this is done, the basic equations become 

( < M)ff + f §e ■*-*(§(+} |f- + 2 f +f 


l-N 


N <P +(?> V 3 ^ + rj 9© q + f 9rj = 0 

N ' » ( ' 3n n 3A n n uj u 


3f 


3rj q 30 q nw 30 


2 1 n N f + (<M) + -q-|e i^=0 


(18) 

(19) 

( 20 ) 
(21) 
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The parameter t\j which appears here is, for the moment, unspecified. 

Next, the boundary conditions must be examined, to see if they are 
compatible with the similarity assumption. At the shock ( fj=l , - 4g- 6. Q ^ + -Ijr) 
defining e Q = 0 , Eqs. (13), (14a), (15), and (14) become 

VQ,e)\[- <PU>e)\ = i 
f(i,e)= 4 >0,e)+ -yjr- 

Co K s 

The first three of these are independent of the time if the initial pressure in 

* 2 . 

the undisturbed medium, '■j} , is small compared with ^ , which is 

of the order of the pressure being generated at the shock. This condition will 
certainly be met whenever the fluid-mechanical model is appropriate. Thus 
the question of whether a similarity solution is compatible with the boundary 
conditions depends solely on whether the form of the internal energy function 
F is such as to permit the time dependence to be eliminated from Eq. (25). 
Sedov^ has pointed out that this can be done wherever the internal energy is 
of the form 

e = <f>(p) (26) 

where y? is any function of the density. For such a case, Eq. (25) becomes 


l*> 


'' 0, e) 


( 22 ) 


(23) 


(24) 


(25) 
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f(!,e) =f>. f 0,6) ip |/3 p 0, ej\ 

( 27 ) 

and all explicit time dependence is eliminated. Thus, a self-similar solution 
is possible whenever the medium obeys the equation of State (26). In this case, 
the boundary values at the shock can be conveniently found by solving Eqs. 

(22) - (24) for <p(j>6), , and in terms of (/,&) 

<P0,e) « fO.e) = i- ( 28 ) 

= T ro7e>) ( 2 9 ) 

When these relations are substituted into Eq. (27), the result is a single 
expression for the density ratio at the shock 

r(ho) = (30) 


Having found ^ (l } d) the other quantities at the shock are found from Eqs. 
(28) and (29). 

5. Perfect-gas Approximation. The target materials of greatest in- 
terest in the present investigation are metals, for which the equation of State 
in the range of pressures appropriate here is the Mie-Gruneisen relation^^ 

e 6 »>p) - (31) 

where the subscript c denotes the cohesive contribution, and where P 
is the Gruneisen constant, which depends weakly on f> . The cohesive 
contributions can be found from measured shock-wave data: along the Hugoniot, 
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Eq. (31) takes the form 


*.w - •.« ■ M '‘XV 


Subtracting this from Eq. (31) then gives 

,„V_ *■- ±a!£L 

" f) pr (p) 

The Mie-Gruneisen equation can be rearranged as 

A 


e = 




- A(f) 


where 


(32) 


Mf) = 


1 

f r (r) 





e n(e) 


(33) 


All of the analyses of this report use only the leading term of Eq. (32), which 
can be accommodated in a self-similar solution. The present solution will 
therefore be valid only when the pressure is sufficiently high that can be 

neglected in comparison with the leading term. In actual fact, most impacts 
will span a time interval during which this approximation fails. It is impor- 
tant to realize that the pressures at which A(/°) is too large to be neglected 
are nevertheless sufficiently high that the compressible-fluid approximation is 
well justified. 

The analyses presented here have made the further approximation of 
neglecting the variation of the Gruneisen factor P (f>) , replacing it by a 

constant, denoted by /- / . Under this approximation, the equation of State 

becomes that of a perfect gas of constant specific -heat ratio / , namely 


(?' Of 


(34) 
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The use of this state equation amounts to a high-pressure approximation of 
the Mie-Gruneisen relation, with taken as a constant, equal to /-/ 

and it makes available all the results of the extensive literature dealing with 
blast waves in a perfect gas. Because P (p) is typically about 2. 0 to 3. 0, 
'Y will be on the order of 3. 0 or greater. It should be borne in mind, 
however, that the similarity solution is not limited to the predictions made 
with the perfect-gas model; the variation of P could be accounted for, but 
is neglected here as a matter of convenience. 

When the perfect-gas approximation is made, the energy equation 
( (19) or (11) ) becomes 


+ u 

dt 




- 0 


(35) 


In terms of the similar functions, this is 


-2 



y£ dy\ u) (9£ 
<// / q \96 


y£ ££_) 

V zej 


^0 


(36) 


It is interesting to examine the accuracy of the perfect-gas model for 
a particular material. The equation of State enters the analysis by way of 
Eq. (11), which states that the entropy of a given particle does not change, 
after being processed by the shock. In general, the entropy of a medium is 
a function of any two thermodynamic variables, for example the pressure and 
density. In making the ideal- gas approximation, the entropy is assumed to 
depend on these only in the special combination 

s = fen (p/p r ) 

To assess the quantitative validity of this assumption, then, it is necessary to 
determine how accurately the isentropic states of a given material are 
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approximated by the relation 


Jt* v-fv — ■/ £# 1 / p - constant (37) 

An explicit equation of state for iron, developed by a group at the Los Alamos 

IS) 

Scientific Laboratory, has been published by Bjork. ' ' This equation has been 
used to calculate the pressure as a function of internal energy at various den- 
sities. The resulting map of thermodynamic states is shown in Fig. 1. Also 
shown is the Hugoniot, i. e. , the locus of states which can exist behind normal 
shock waves. This locus is the set of points which satisfy the condition 


e - 


A_ 




(38) 


where p o is the normal density of iron, 7. 86 grams per cubic centimeter. 

The experimental Hugoniot curve, recently determined up to pressures of nine 

(g\ 

megabars by Altshuler et al v ' is included, and agrees quite well with the the- 
oretical curve. In addition to this experimental check, the pressures shown on 
Fig. 1 at zero energy' agree quite well with the Mie-Gruneisen values 


* 


iA = p r(p) A(p) 


(39) 


given by Altshuler et al in an earlier publication. Thus Fig. 1 may be con- 

sidered a valid representation of the thermodynamic states of compressed iron. 

Three isentropes have been added to this map, by a trapezoidal- rule 
integration of 


or 


'P/Po 

(e/?oY 


(40) 


It should be noted that Eq. (39) provides a simple way of calculating the relative 
contributions from the two terms in Eq. (32) 


e 


^ r pAjp)r'(p) 

pF(f>) L ^ 


pP(p) ^ 
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which holds for constant entropy. The three isentropes shown may be 
thought of as describing the histories of three particles which lie along 
the axis of symmetry in a given impact, and which are therefore processed 
by successively weaker stages of the same shockwave. Particle 1, for 
example, is raised from its normal state to the point shown on the Hugoniot 
by a shock wave moving at 20. 8 km/ sec; its subsequent expansion to low 
density takes place along the isentrope shown. In similar fashion, par- 
ticles 2 and 3 are processed by the shock when it is traveling at 16. 5 and 
8.62 km/sec, respectively. The -p , p -coordinates of these three isen- 
tropes are re-plotted in Fig. 2. Large portions of these curves are accurately 
reproduced by a constant value of / , i. e., they have a constant slope. 

6. Procedure for finding y* , In applying a constant - theory to 
any given impact problem, only a single value of Y may be used. If, in 
the example cited above, one chooses the value appropriate to particle 1 at 
the shock, the approximation will begin to deteriorate for the lower-density 
states of this particle, and will be less accurate for all states of particles 
2 and 3. In order to minimize the error, it would appear that Y should 
match the higher-density, higher-pressure portions of the flow. Partly for 
this reason, and partly for ease of application, the procedure recommended 
here is to choose Y so as to match conditions at the impact point. For 
any impact, one may imagine that there is a small region in the immediate 
vicinity of the impact point, where the collision is equivalent to the planar 
impact of two semi -infinite bodies. 
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We shall see below how to calculate exactly the pressure ■p'f , internal 
energy , density p t , and particle velocity U , t , behind the shock 

wave that starts into the target, as well as the speed of the shock, ot s . 
For a strong shock in an ideal gas, use of the state Equation (34) in the 
strong-shock form of the Rankine-Hugoniot relations, Eqs. (27) - (30), 
reveals that these quantities would be related by 


- ±L „ _ Z z 

U s - - rn ’ Po - r-i 7 • (7+1)' 5 


(41) 


The procedure suggested here is to choose the value of Y that satisfies 
Eq. (41), wherein the quantities , e ; , ~p' l , and p t used are the 
actual values that would occur at the impact point in the target material. 


For example, if the target is compressed to twice normal density at the 
impact point, the "effective ■/ " for the process would be 3. 0. 


It should be noted that ? can be found from any of Eqs. (41), i. e. , 


7- Z 


-I = 2 


a u t . _ * 1 


p.fa - 1 


Zu A 

~ '2A — 1 


Fig. 3 shows the effective Y for iron, as a function of shock speed U s . 


This curve was taken from the calculated Hugoniot of Fig. 1, together with 
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Typically 


the experimental data of Altshuler et al^ and of Walsh et al. 
of most materials, the effective Y is on the ordei of 10 to 20 for weak 
shock waves, and is more like 2 to 3 for strong shock waves. 

To calculate the actual conditions at the impact point, only the Hugoniot 
curves of the projectile and target are needed. Consider the impact, at speed 


Immediately after impact, a shock wave moves into the target at speed , 

and the material behind it, at pressure -fl>£ , follows the shock at speed Vf 




V , of two semi- infinite bodies 



Meanwhile, another shock wave propagates back into the oncoming projectile 

A 

material, reducing its speed from V to V. , and raising its pressure to 

-p/j . Equality of pressure and particle speed across the interface requires 
A A 

that -ja>. - and V • = \/. . Consider now the shock propagation 

in each medium separately, in each case using a coordinate system such that 
the shock moves into a medium at rest 
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The problem is to choose, for a given impact speed y , a pair of values 

A A a 

V and v. such that p. - P. and V- = y. . Olshaker and 

S ' C t C t t 

Bjork^^^ have described a convenient graphical method of solution. For 
the present purpose, where the solution is required over a range of impact 
speeds, it is much simpler to turn the problem around, choosing the pressure 
at the interface, and asking what impact speed this corresponds to. The 
choice of p. (which equals p. ) determines the particle velocities U f 
and , from the Hugoniot curves for the two materials. The impact 

speed is then determined by the equality of particle speeds in the laboratory 
frame of reference, namely 

V. = U. - v. = V - U t 

t 1 v i 1 


A (43) 

or V - + u t 

The value of Y can then be found from any of Eqs. (42). This process 
is illustrated in Fig. 4; here are shown the -p l/ , u t Hugoniot data for 
iron and lead, taken from Altshuler et al. ^ Choosing an interface pressure 
of 8 megabars gives = 6. 6 km/ sec and 7. 0 km/ sec in lead and iron, re- 
spectively. Thus the impact speed required to generate these conditions must 
have been 13. 6 km/ sec. If lead is the target, the effective Y is determined 
as follows: 

a) The shock speed is determined, either from a graph of U. t ys, u s , 

or vs. U 5 , or else from Eq. (14a) 

U «. = 3^2 - = 10. 69 km/ sec 
5 fo“-> 

b) Y is then found from Eqs. (4) 

^ _ / = 2 . 24 

U, 

Repeating this process at a succession of values of p determines tfiv) , 
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as shown in Fig. 5. 

The determination of -tf (]/) takes on a simple form when the shock 
velocity and particle velocity are linearly related 


U 5 - C -h su, (44) 

The measured Hugoniots for many materials are well approximated by this 
relation, at least up to pressures around two megabars. Table I lists values 
of C and S for a number of materials. Using Eq. (14a), the particle 
velocity corresponding to a given pressure is 


U, - C 


//Hs ft/ftC 2 -/ 
25 


(45) 


The impact speed can then be found by applying this formula to the projectile 
and target materials, and Y is found in the usual way, for example from 

^ ^ ^ ^ u>i) I ( 46 ) 

The linear relation between shock velocity and particle velocity ceases to be 
valid at extremely high pressures, but our present experience indicates that 
it is a satisfactory extrapolation method for impact speeds up to 50 km/ sec. 
The Gruneisen constant, f 1 , may be as large as 3. 0 for the materials of 
interest here. Thus, as long as Y is less than 4. 0, the perfect-gas approx- 
imation may be interpreted in the sense indicated above, namely that the 
leading term of the Mie-Griineisen equation dominates, and the choice of y 
is essentially an approximation to the Gruneisen factor. If the "effective 
Y M is greater than 4, however, no such interpretation can be made, and the 
use of a perfect gas must be viewed as an attempt to match the entire Mie- 
Gruneisen equation with a single term. Obviously, the approximation will 
not be as good in this regime, and the theoretical predictions at low impact 
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speed (i. e. , high ) must be considered less reliable on this account. An 
estimate of the impact conditions for which Y is less than 4. 0 can easily 
be derived from the relations above, from which it can be shown that Y will 
be less than 4. 0 whenever the impact pressure is greater than 

A > po c ^ ^ 

As s ranges from 1 to 1. 5, this pressure varies from f> o 0 to 

Thus, the present theory may be expected to agree with experiment at lower 

£ 

velocities, the lower the value of p C . Reference to Table I, for example, 
shows that, at a given velocity, impacts into lead would be better described 
than those into, say, copper. We shall see below that this is indeed the case. 

7. Conservation of Energy and Momentum. The total energy and 
momentum of the system must be conserved, as may be confirmed by forming 
the proper volume integral of the vector equations of motion, Eqs. (1) - (3). 

The actual integrals, whose values must be constant, may be derived as follows 
consider, as the mass element, a ring of volume r ci fd 0 * 2 ft' r Jovo 0 



The total energy £ and momentum P are 


E = 


e 


i 2 2 

\d i-Axr 


j J p ' 2irr Ain-' 6 d r U 0 


(47) 


= Zfr 'po A £5 f ^ & f ' (f + J ]) t 1 %aL 7 ^ 6 
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( 48 ) 


n p (£} 

(u -surA*i Q) p ■ 2 ?rr 2 Q dr e/G 


. 3 nv pn.(©) /. \- 

= Z ft f> 0 ^UyfV 6 (<j> 'CO& 6 - CD MwG ) p rj Z ffC fj ei G 

X 

Here we encounter a fundamental difficulty. If we are to have a self-similar 
solution, the differential equations require R g - A i . However, a single 
value of N will not permit both of the relations above to be independent of 
time. Constancy of energy can be achieved only with N =2/5, while momen- 
tum conservation requires J\j = 1/4, and in either case the parameter A is 


used to match the quantity being conserved. Thus it appears at first glance 
that a satisfactory solution cannot be achieved under the assumption of sim- 
ilarity. The essence of the difficulty is that, having used M to make one of 
the integrals independent of time, only a single free parameter, A , is left. 
But we still have two quantities to be matched, as well as a second integral to 
be made independent of time. We will describe, in Section C of this part, one 
method for overcoming these difficulties. The essence of the method is that 
N is determined by a totally different consideration, and a second free 
parameter is introduced in such a way that both conservation conditions may 
be satisfied simultaneously. 

8. Final Form of the Problem Posed. The final set of differential 
equations that must be solved consists of Eqs. (18), (19), (20), and (36), sub- 
ject to the boundary conditions at the shock r| s I , — • Tl jZ -06+ ft/ Z 

4 (!,$)• f 0,e) ‘ -/j- > <>(I,9)~0 ; 40,e)=~ < 49 ) 

The value of Y is to be found, for each specific impact case, from the 
method outlined in Section A-6. The specification of N will be described in 
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detail below. 


As a final note on the mathematical nature of the problem, it should 
be pointed out that two pairs of characteristics of Eqs. (18), (19), (20), and 
(36) can be found by standard methods. One pair corresponds to the particle 
paths of the original unsteady flow, while the other pair, related to the Mach 
lines of the original flow, reveals that the equations have elliptic or hyperbolic 
character, according to whether 

i 0 (50) 

The line in the Q plane along which this quantity vanishes is referred to 
below as the "sonic" line, because of its relation to the Mach lines of the 
original flow. 
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B . Spherically Symmetric (Constant-Energy) Solution 

This section describes a solution in which the shock propagation 
into the ta,rget is represented by one half of a spherically- symmetric 
disturbance: I 



Such a solution allows variations only in the radial direction, and does not 
describe the pattern of mass ejection from the expanding crater. Moreover, 
we shall see below that only the .total energy may be conserved, but not the 
total momentum. The justification for the use of such an apparently defective 
model is that it closely approximates the results found from a more accept- 
able (asymmetric) model. In this section, we present the symmetric solution 
without apology, deferring its justification until after the asymmetric solution 
has been treated in Section C. 

1. Solution for R s (t) . When the flow is spherically symmetric, 
all derivatives with respect to Q vanish, and the similarity equations be- 
come ordinary differential equations. Denoting the ordinary derivative with 
respect to rj by a prime, these are 

y>' * + zjf) = o 

-■4r * + (*-<))*' + yr = o 

- zJ irf *')-o 


(51) 


(52) 


(53) 
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These maybe solved explicitly for the derivatives , <fi' , and -p 'in 
the form 


^ + Z bJ±f _ 


1 


M 


N 




0-1) 

^ *0-1)+ 

|0‘l) y J 

I L l£-Z 

V 0 

/-n f 
N y 

_ 4 

(*-i) 

0-r 

fzizi 
L n 

) 2 - ^ 

; v* 

n J 

- * 7 

\ 


( 54 ) 


(55) 


rf 


(M) - 7 ? 

The boundary conditions at the shock are 


(56) 


*(')= PO). ; <f>(0 


(57) 


Equations (54) - (56) (with ^ = 2/5) were first presented by G. I. Taylor, 
who worked out a few numerical and approximate analytic solutions, for y 
ranging from 1 . 2 to 1.67, the range appropriate for gases. Subsequently, 
an analytic solution of these equations (also with M =2/5) was published by 


(13) 


J. L. Taylor/ 14 ^ Latter/ 15 ^and Sakurai/ 1 ^ Simultaneously with Taylor's 


( 6 ), 


work, Sedov had also found this analytic solution. He worked in a different 
coordinate system, where 


V 

1 


^ _ N 

f -~ 


rf 


y f 

In these coordinates, the basic equations may be combined to yield a single 
differential equation involving only ^ and (see Appendix A for details of 


(58) 


the derivation): 
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i-N d. 


& =-(r-0 + 



]* 

■-Nf-y] 

IL 

-(5 r-i) 


r 

o- 

W)] 

(rt-N)fa-l) 


7 d* 


Having y, as a function of , the quantities fj and p are found by 
integrating 

^ (# -A/) l - 


1 * 




The boundary conditions at the shock are 

ZH mm _ 2rd % (r-i) 


- 

SH (r+0 


SH 


(?+>)' 


The parameter SI must be specified, before solutions of these 

equations can be found. It appears that physically acceptable solutions of 

these equations exist only when SI = 2/5, a value which conserves the total 

energy, according to Section A- 7 above. When SI is taken to be different 

from 2/5, the solution exhibits infinite slopes. In fact, by working in the 

coordinates, the solution can be shown to become double- valued. 

Figure 6 shows results typical of those found in the range 0.25 £ SI — 0.4 

when a solution of this sort is attempted. This nonexistence of symmetric 

solutions apparently explains the difficulty encountered by Davids, Huang, 

17 

and Juanzemis in attempting to find a spherically- symmetric solution for 
constant momentum ( Si = 1 /4). 

In what follows, SI is chosen as 2/5, and the terms "constant-energy' 
and "spherically- symmetric" are used interchangeably in referring to the 
solution. 


( 59 ) 


( 60 ) 


( 61 ) 


( 62 ) 
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The closed-form solution found by Sedov (or equivalently, with q as 
independent variable, by Latter, J. L. Taylor, and Sakurai) is 


, , r(r-l) ,z z/5 - * 

'%(*) - — 2 — * 


- Z/S 


and expressions for the remaining details of the solution can be found, for 


example in Sedov^or in Hayes and Probstein.^^ Detailed numerical results 


are given only for Y in the range appropriate to gases. Calculations have 
been done for Y =2,3,4,5,6,8,10,12,16, and 20, and the distributions of 
velocity, density, and pressure are shown in Figs. 7a-j. These figures dis- 
play the usual feature, that in all cases the density drops off rather sharply 
behind the shock, indicating that most of the mass processed is 
concentrated near the shock. As Y approaches seven, the distributions 
approach the exact solution, for Y - 


> 


<p = ir n 




For Y greater than seven, a cavity begins to form at small values of 

( 6 ) 


j' 8 * (e ♦ i • 2wrVr * p fci J - + t “ ! ') 2 v rVi 


= Z */°oRs R s 


( 63 ) 


(64] 


as pointed out in Sedov, and the particle velocities show a marked increase 
near the edge of this cavity. 

Having these distributions, an explicit description of the shock prop- 
agation can now be given, if the total energy of the system is specified. 
The sum of the internal and kinetic energy of the fluid that has been set into 
motion is given by the integral 


(65 
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where, 


x > w ■ {' (pby t (66) 

This integral has been evaluated for the values of Y mentioned above, and 

Cl 3) 

is shown in Fig. 8. The values obtained by G. I. Taylor at lower values 
of Y are also shown here, together with the result reported by Davids and 

Cl 9) 

Huang 1 'at Y = 16, about which more will be said below. A further check 
is the value X ( =1/72 when Y = 7, which can be found by a simple 
integration. 

If the total energy £ is now specified, a simple differential equation 
for R s (t) results 

E = 2 irp 0 R* s R* I ,(/) (67) 

The term is three times the target mass processed up to the time t . 

Thus, 31 (?) may be thought of as a dimensionless coefficient giving the ratio 

# 2 

of the mass-averaged value of e+^-to the quantity £ , i.e. 



• z 

But J? is proportional to the energy at the shock 


^k±<l 




Thus 




4 L»» 

3^/) z (e 

\ / SH 


(69) 


(70) 
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Because most of the mass is concentrated near the shock, the mass- 
averaged value of any quantity is very nearly its value at the shock. Thus 
the factor is a good approximation to X, (?) , as shown in Fig. 8. 

This factor originates from Eq. (69), which states that the larger the value of Y , 
the larger must be the shock speed, if a given energy is to be achieved behind 
the shock. We may attach the same significance to I | (V) : if a given energy 
is to be distributed in two materials for which the /S differ, the shock speed 
will have to be greater in the material having the larger Y . 

The solution of Eq. (67) is simply the Taylor solution for a strong blast 


wave 


K.tt) = 


Z5 





(71) 


Here the influence of / is shown more clearly. For given E and , 

the shock radius will grow more rapidly for large values of : 



To apply Eq. (71) to a given case, the total energy E must be specified. 
In all the applications made below, this energy is taken to be the kinetic 
energy of the impacting particle. 

Figures 9ct and 9b present a comparison of the Taylor solution with 
experiment. Eichelberger and Gehringf^ have published time histories of 
the shock and crater radius, for a Lucite block struck by an iron pellet at 
4. 6 km/sec. These time histories are compared with the prediction of 
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Eq. (71) in Fig. 9a. The agreement is reasonable, although the measured 
rate of advance of the shock is faster than that of the constant- energy solution, 
corresponding to a value of f\| around 0.67. It is important, in this connec- 
tion, to recall that any constant-energy, self- similar solution will have 
N =0.4. Changing the value of "/ will change only the amplitude. In 
addition, any equation of state that permits a similarity solution will also 
leave M unchanged in a constant-energy solution. Furthermore, in a solu- 
tion which conserves both energy and momentum, the evidence advanced in 
Section C below suggests that, if anything, will be less than 0.4. On the 

other hand, it is true in general that the shock speed will decelerate from 
the blast-wave behavior (j? 5 «—£^)at early time to the acoustic limit t) 

at late time. The data of Fig. 9a apparently lie in a transitional regime 
between these two limits. Thus it appears that the only way to achieve a 
better comparison with the data of Fig. 9a is to include the nonsimilar effect, 
associated with the term in the equation of state. The present report 

does not treat this effect, but it is encouraging to note how well the blast- 
wave theory does in spite of this deficiency. Moreover, the blast-wave 
theory can be expected to improve at higher impact speeds. 

Figure 9b presents a different type of experimental check on the blast- 

wave theory. Here we show the wave speed versus shock position, measured 

(2) 

in a very ingenious series of experiments by Frasier and Karpov. ' The data 
pertain to a wax target, struck by an Ethocell pellet at 4 km/sec, and 
apparently lie in a range where the shock strength is too low to justify the 
strong- shock assumptions. However, the data do not appear to be incon- 
sistent with a transition from a high shock speed down to the stress-wave 
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velocity. The constant-energy solution is not extrapolated beyond this 
velocity, since the blast-wave approximation is clearly inadequate in that 
range . 

In general, then, the blast-wave solution may be considered a rea- 
sonable approximation to this very limited amount of experiment, and, if 
used with caution, may be expected to serve as a suitable basis for crater 
prediction. 

2. Crater formation criterion . The blast-wave solution given in 
Eq. (71) may be considered a valid description of the shock propagation in 
the target so long as the pressures being generated are large enough to 
justify an inviscid- fluid model. These pressures decay rapidly, however, 
and to provide a valid solution at later time, a transition is needed from the 
blast-wave model to one which properly describes the plastic flow, and 
ultimately the elastic response of the target. Such a solution would predict 
the configuration in which the target material finally comes to rest, and 
the definition of the final crater dimensions wd uld be unequivocal. A solu- 
tion of this sort is not presently available, however. In its absence, the 
best that can be done is to identify the point at which the transition from 
blast-wave theory ought to occur, and to make an estimate, based on con- 
ditions at that instant, of what the final crater dimensions will be. 

The method for predicting crater size that is adopted in the present 
report is that the shock radius at the instant at which the blast-wave solution 
is to be cut off is the radius of the crater that will ultimately develop. The 
use of such a procedure is equivalent to the statement that all of the material 
processed by the high-pressure phases of the shock wave will ultimately be 
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ejected from the target. 

The cutoff point must be the point at which the shock-wave intensity 
has decayed to some preassigned level. If we require that the pressure 
behind the shock shall have fallen to the level fP (to be identified below 
with the material strength), then, since the pressure behind the shock is 
given in general by 


*/ = 77T &■ *• 

we would identify the crater radius as the value of 


rj- / p 


2 

If the solution for the shock radius is written as 

c \i 5 


f? s = A t ^ 5 > 


A = 


25 




(72) 

R s at the instant when 

(73) 

(74) 


then the time can be replaced in favor of the shock radius, according to 

* = (%S ( 75 ) 


The shock speed is therefore 




3/a 


(76) 


Equating this to Eq. (73), and replacing J? s by the crater radius , we 

have 


2 

5 


a /A_f - JUTm. 
A U. / ' / 2 ft 


or 




Z = 


c \Tr(r+\)l(r) F, 


(77) 
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Taking £ to be the kinetic energy of a spherical projectile of density 
and diameter ci. then gives 


This solution is shown in Fig. 10 for various values of y . For a given 
projectile -target combination, y is large at low impact speeds, and de- 
creases as V increases. Thus the crater radius is predicted to grow with 
a power of Y somewhat less than 2/3. In most cases, the power is approx- 
imately 1/3. The influence of y displayed in Fig. 10 is the same as that 
mentioned earlier: at a given kinetic energy, the higher values of y produce 

faster shock waves, which will penetrate more deeply before they decay to the 
prescribed pressure IP . Stated another way: if a given energy is to be 

added to the target by a shock wave across which the energy change per unit 
shock speed is relatively small, then the shock speed itself must be relatively 
large. We will return to a discussion of this point below. 

The value of IP must be chosen, in order to predict actual crater 
dimenions. It was pointed out above that the inviscid-fluid approximations 
are valid only so long as the pressures being generated at the shock are 
large compared with the target's resistance to shear deformation. Thus, 
the pressure at which the blast-wave solution is to be cut off must be a 
measure of the shear strength of the material. Part II, below, points out 
that the proper value to use for this property is the. intrinsic shear strength, 
which lies between the limits Q,ji0 and QrjT.lt' > G- being the dynamic shear 
modulus, as measured, for example, by ultrasonic techniques. Using 
these values for the intrinsic shear strength, the predictions of Eq. (71) 


^ _ / ! fd 

ob (/2 (rrl) !,(>■) P 



RM-1655-M-4 


32 


are compared with experiment in Figs. 11-13. In preparing the theoreti- 
cal curve for iron striking lead (Fig. 11), the variation of / and )/ 
was taken from Fig. 5, and the material properties used were found from 
Table 4 * . The lower strength leads to a deeper crater prediction, since 

the shock wave must penetrate more deeply in decaying to a lower pressure. 
The experimental data, taken from the compilation by Herrmann and Jones^^ 
are seen to agree quite well with the higher- strength prediction. Figure 12 
gives results for iron striking iron. In addition to low-speed experimental 
data, this figure also shows Bjork’s machine solution. Again the agreement 
is reasonable. Figure 13 shows the results found for aluminum striking 
copper. Here the agreement with experiment is somewhat poorer. The 
Hugoniot data for iron, lead,and copper were taken from Altshuler et al,^ 

while for aluminum the low-pressure data of Walsh et af ^ were used, together 

121 ) 

with the high-pressure estimate made by Lake and Todd. 

On the basis of this evidence, it would appear that the best agreement 
with experiment is achieved by using the value GrfZff for the intrinsic shear 
strength. Part II, below, presents a discussion of why such a conclusion 
might be expected. 

The approximations on which this theory is based will improve at 
higher impact speeds. Thus it is not surprising to find such good agreement 
for the case of lead targets, in which hypervelocity conditions are achieved 
at relatively low speed. The fact that the slope of these data and of Bjork's 
results are well matched lends substance to the belief that crater radius will 
grow somewhat more slowly with impact speed in the very high-speed range, 
compared with its growth rate at lower speed. 
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3. Scaling Laws. The only target properties that appear in the 
crater- size formula are the strength JP and the equation of state, as 
represented by X^/)- For the materials considered, however, the latter 
influence is of secondary importance, since y(V) is approximately the same 
for a variety of projectile-target combinations. This point is illustrated in 
Figs. 14a, b, and c, which show the functions y(V) for targets of Aluminum, 
Beryllium, Iron, Molybdenum, Columbium (Niobium), Vanadium, Tantalum, 
Tungsten, and Lead, being struck by projectiles of Fused Quartz, Aluminum, 
and Iron. These calculations used the data of Table I and the equations of 
Section A6, Part I. Crater- size predictions for these projectile-target com- 
binations can now be made by using these results for /(V) in conjunction with 
Fig. 10. Thus an important scaling law is that the crater radii produced in 
two targets, by a given projectile impacting at a given speed, will have the 
ratio 

M. . 


jp. 


. 1/3 


Sol 


p. 


(79) 


The fact that 7 ^(V)is not greatly different for various materials also permits 
a simple scaling with respect to projectile density. For impact into a given 
target at given speed, the craters made by different projectiles are related by 
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/ x \'J» 
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(80) 


Finally, it should be noted that the scaling shown here differs from 

( 12 ) 


that recommended by Olshaker and Bjork'"^ Their scaling is based on the 
observation that craters produced in a given target by projectiles of various 
materials, all having the same mass and velocity (and thus different diameter 
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are proportional to the particle velocity generated in the target at the impact 
point. Thus, for B striking A , 


'U 


= con sit x 




(81) 


B— '/ in /\ , at the impact point 
The experimental data which provide the basis for this scaling are all taken 
at relatively low impact speeds. Because the particle speed at the impact 
point is one half the impact speed for like-on-like impacts, the scaling law 
can be written as 


& 


u >). 


B— *A 


R«) 


A->a 


v/z 


(82) 


If the shock driven into A by B is stronger than that driven into A by A , 
this scaling would predict a larger crater for B striking A than for A 
striking A - The present theory gives a different scaling law; for given 
kinetic energy of the projectile, and given JP , 

V/3 

(83) 






A — > A 






8— 


It can be seen from Eq. (46) that Y decreases as u t increases. Thus, 
for the case mentioned above, where U ^ > S/Jz > it follows that 


7l < / 


8->A 

, so the prediction of the present theory would be 1^) <|? \ 

* '*->A C A— 


just the reverse of that predicted by the Olshaker-Bjork formula. The reason 
for the difference is simply due to the influence of y that has been discussed 
above. Although the present theory uses the shock strength at the impact point 
to determine •/ , nevertheless the solution does not match the shock speed at 
that point. Rather, the role of Y is to characterize the strength of the shock 
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throughout its entire subsequent history. The stronger the shock, the 
less rapidly it must travel, on the average, in order to transfer a given 
energy to the target. 

Unfortunately, the data required to resolve this difference are not 
presently available. In actual fact, of course, the speed of the shock as 
it advances into the target will remain close to the one-dimensional, impact- 
point value until the time when the projectile has been destroyed. Thereafter, 
its speed will be governed by the requirement that the total energy and momen- 
tum be conserved. One must certainly expect that the Olshaker-Bjork type 
of scaling will be correct when the entire cratering process is controlled by 
conditions during the time when the projectile is being destroyed. This 
condition has been shown to exist at low impact speed. Whether it will con- 
tinue to hold at higher impact speeds is presently open to question. On the 
other hand, the blast-wave theory in its present form is valid only at impact 
speeds such that the collision time is a small fraction of the time required 
for the crater to form, and so its predictions are relatively insensitive to the 
details of the impact process. The fact that it leads to a scaling law differ- 
ent from that observed at low speed is therefore not surprising, but must 
be regarded as tentative, pending the availability of a less approximate 
theory, and of definitive experiments. 

It would appear that these qualitative conclusions would remain un- 
changed even if more realistic Gruneisen coefficients were employed in a 
blast-wave solution. The extent by which they might be modified by the re- 
sults of a non- similar solution will form an interesting area for further 
research. 
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4. Comparison with Other Theories . It is of interest to compare 
the present analysis w'th a similar theory given by Davids et al.^^’^^ 

These authors also use a spherically symmetric model, and represent the 
equation of state by a gas of constant Y , which they determine by fitting 
the Hugoniot (not the isentropes) by an equation of the form t^? p . They 
work out a numerical solution for Y =16, and N =2/5, but apparently 
do not take advantage of the closed-form solution. The value which they 
find numerically for the cavity radius differs slightly from that given by 
the exact result, and this small discrepancy may explain why their value 
of !,(/,), plotted in Fig. 8, is somewhat below the present results. As noted 
earlier, they have also attempted a solution for constant momentum ( N = 1/4). 
The principal difference between their theory and the one described here lies 
in the criterion used for crater formation. Like Bjork, these authors do not 
incorporate the material strength. Rather, they define the crater as the cavity 
radius at the point where a certain graph of shock radius versus time is judged 
to have reached an essentially constant value. They infer from this that 
crater radius will vary as the 2/5 power of the impact speed. 

The present formula for crater size is very similar to that derived 
empirically by Eichelberger and Gehring.^ Their result states that the crater 
volume varies directly with the kinetic energy of the projectile, and inversely 
with the Brinell hardness of the target, B 

-§- #•*>/ = -§- 

where E and are taken in cgS units, and B in the customary units of 
kilograms force per square millimeter. The present result is 
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Z 

3 



3 (7+1)1 t (f) 


£ 

G/Ztt 


which differs from Eichelberger and Gehring's formula by the influence 
of y , and by the use of G/Zlf . There appears to be no simple relation 
between the dynamic shear modulus and the Brinell hardness, so that no 
general comparison can be made. For the case of copper, Gr/zif is about 
ten times the Brinell hardness, and the two crater-size predictions are equal 
at y = 4.6, a representative value. In general, both predictions yield the 
same order of magnitude, and differ principally in that the present one in- 
dicates a lower growth rate of crater size with impact speed. 
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C. Asymmetric Solutions^ 

In this section we present a blast-wave solution which explicitly des- 
cribes the pattern of mass ejection from the target. Because of the lack of 
spherical symmetry, two independent variables enter, and the solution of 
the problem is considerably more difficult to find. Certain approximations 
are resorted to, in order to obtain partial solutions. From these, a very 
important result is found, namely that this more difficult, but more physically 
realistic solution is for practical purposes the same as the vastly- simpler 
constant-energy solution described above. 

In addition, this section returns to the question of the simultaneous 
conservation of energy and momentum, and describes one method by which 
this can be achieved. The method is illustrated by applying it to the case of 
a one-dimensional impact, such as might occur when a thin disk strikes a 
target. 

The first six of the subsections below describe the axisymmetric case, 
while the last three present the one-dimensional solution. 

1 . Two-Dimensional Solution. The model which allows for spatial 
variations in two directions behind a hemispherical shock wave was intro- 
duced in Section A, above. The similarity equations, Eqs. (18) - (20), and 
Eq. (36), are partial differential equations, containing both f| and Q 

❖ 

The term "asymmetric" as used here refers to symmetry with respect to 
the target surface, while the term "axisymmetric" refers to symmetry about 
the axis along which the projectile impacts. All of the solutions of this 
report are taken to be axisymmetric, but may be either symmetric or asym- 
metric with respect to the target surface. 
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as independent variables. They are of mixed character, containing both 
elliptic and hyperbolic regions, and, furthermore, they must meet a zero- 
pressure boundary condition along a line whose location is unknown in 
advance. To make matters worse, the differential equations contain a para- 
meter /\] whose value is unspecified. No attempt has been made to solve 
these equations. Instead, partial solutions are sought by restricting 
attention to conditions along the axis of symmetry. In this way, we can 
learn a great deal about the solution, with relatively little effort. The 
most important item uncovered is the criterion for choosing M : for each 
y , M must be chosen so as to permit a smooth transition from the el- 
liptic to the hyperbolic region. 

In the similarity coordinates, the flow field has the appearance 


In the undisturbed region, it is as though fluid particles were all converging 
radially toward the origin. They pass through the shock, and are ultimately 
ejected toward the low-pressure region outside. On the basis of the spherically 
symmetric solution, we may expect that most of the mass processed at any 
instant is heavily concentrated near the shock. Thus the problem bears a 
marked resemblance to the steady, hypersonic flow over a blunt body 



L 
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I 

and we may expect that analytical methods which are successful in treating 
that problem may also be used to good advantage in the present case. One 
such method is to investigate the solution along the axis of symmetry: 

9=0: 0 ^ — / 5 0 = TT , O rj ± oo (84) 

Along this line, the axial symmetry of the problem requires that all first 
derivatives with respect to Q be zero, except for the derivative of the tan- 
gential velocity component. This component is antisymmetric in 0 ; 

for example, in the case where material Is flowing outward from the region 
near the axis 

denotes resultant 
velocity: 

• ^ = U + xtr 

the distribution of &) at a given radius has the appearance 





Thus — — (rj ,0)^-0 • Denoting this quantity by f : 

r (H) = 2 |f- (f| ; 9 = 0 or it) 

The similarity equations, along the axis of symmetry, become 

(<f>-r\) l t' l + { l' (<t>‘ + Z-j-) =0 


(85) 


( 86 ) 
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I-N 

N 


(87) 


<t> +(4>- 1 )) <P' + -£-=o 

-2if f ~ljr t') = 0 < 88 > 

where primes indicate ordinary differentiation with respect to fj . These 
equations are referred to in what follows as the "centerline" equations. 
Except for the presence of f in (86), these are identical with the Taylor 
equations for a spherically- symmetric disturbance, discussed in Section B. 
The function T’(rj) represents the influence of off-axis conditions, as must 
be expected whenever a partial differential equation is specialized to a single 
line in the plane of its independent variables. 

Equations (86) - (88) may be solved explicitly for the derivatives in the 



The boundary conditions at the shock are 

MO* t(0= i ^(0 - ~j ! r(l)=0 (92) 

These equations contain two unspecified quantities, N and T'(if)) . The 
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next two subsections discuss the determination of these. 


2. Criterion for Choosing N 


The centerline equations obviously 


z y n 

have a singularity at the point where the denominator (</>-(]) — ■ vanishes. 

This quantity is the special case, for co = 0, of the function discussed in 
Section A 8, whose sign determines whether the partial differential equations 
have elliptic or hyperbolic character. Thus the point on the axis of symmetry 
where this denominator changes sign corresponds to the intersection of the 
"sonic line" with the axis. In order that the solution may pass smoothly 
through this singular point, the numerators of Eqs. (89) - (91) must also vanish 
at this point. A little algebra shows that this condition may be achieved simul- 
taneously in all three numerators if 




I-N 

N 


<P' 


(93) 


1 n rj * 

where the asterisk denotes conditions at the sonic point. The function r fl) 
cannot be chosen arbitrarily; thus the only parameter that can be used to guar- 
antee a smooth crossing of the sonic point is M , and this consideration forms 
the criterion for the choice of N . For each Y , and a specification of 
r , n is chosen so as to provide a continuous transition through the sin- 
gularity. Thus M will in general be a function of ? . It should be noted in 
passing that this problem never came up in the spherically- symmetric, constant- 
energy case. There the vanishing of the denominator always coincides with 
either the origin f| = 0 (for f 6: 7 ) or with the edge of the cavity (for Y>1), 
so the entire flow field is elliptic in a constant-energy solution. 

3. Approximations for r(q) . In order to actually carry out a smooth 
crossing of the sonic point, Eqs. (89) - (91) must be solved for various values 
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of N (and given Y ) until such a crossing is found. Before such an inte- 
gration can be done, however, T(l^) must be specified. In actual fact, no 
rigorous determination of ?Y*]) > and with it N (Y) , can be made with- 

out solving the full partial differential equations. Approximations to N 
may be found, however, by approximating Y , and then integrating Eqs. (89)- 
(91). Rather than approximating T itself, one may instead relate T to 
other physical quantities which may be approximated more easily. In par- 
ticular, by differentiating Eq. (20), the 0 - component of the momentum 
equation, with respect to Q , and by then specializing to the axis of symmetry, 
one finds 




/ 


Y± 

de* 


(f\) o) =0 (94) 


from which it is seen that approximations to the pressure distribution can 
be used to generate corresponding approximations to • This process 

can be continued, of course, by taking higher-order derivatives, with 
respect to Q , of any of the equations of motion. Each of the resulting 
expressions will contain at least one unknown function, so the utility of the 
procedure is dictated by one's ability to approximate the unknown function. 

For this purpose, Eq. (94) is especially useful. At the shock, the pressure 
is uniform, while behind the shock it begins to decrease. The rate of decrease 
is faster near Q - - > as the influence of the vacuum outside the develop- 

ing crater makes itself felt. Qualitatively, the pressure distribution would 


Approximations to the pressure distribution are what make possible center- 
line solutions of the blunt-body problem. 
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be expected to have the appearance 


f i 


/ 



z 

} 

L 

2 fj DECREASING 

g*r , . * 

The quantity which is essentially the curvature of these lines at 


0 — 0 » will be zero at the shock, and will become negative with increas- 
ing magnitude as falls below one. Such considerations suggest the 

approximation 


9 V CL 

-J£T (n>0) = -K(l-n) ,0) (95) 

where K and a. are constants. Crudely, one may think of this approx- 
imation as fitting a cosine variation to the curves above, with a multiplicative 
function of rj introduced in such a way as to guarantee zero curvature at 
the shock. 

The constants K and CL must be chosen so as to yield values of 
T which are at most of unit order. This consideration is derived from the 
fact that, near the axis, the velocity vectors are expected, on physical grounds, 
not to diverge very rapidly from the axis 



The angle which the velocity vector makes with the radius is of order 
g — - (*]) 0 ) - 0 — . But both 0 and (f) are of order one or less. Thus the 

order of dti)/dQ must be the same, if the velocity vectors are to diverge 
from the axis at a moderate rate. 

4. Results for a, = 1, K =1, 10 . A limited number of solutions 
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have been calculated, using the values a. = 1, /< =1.0 and 10. For a 

given value of Y , and selected values of N , Eqs. (89), (90), (91) and 
(94) are integrated by a Runge-Kutta procedure, starting from the shock 
values given by Eq. (92). The results found with a. = 1, K = 10, 7 =3, 

and N| =0*374,0. 375, and 0*376, are shown in Fig. 15a. The first graph on 
this figure displays the denominator, whose vanishing identifies the H sonic n 
point. The second plot shows the function H(j ]) , defined as 


which, according to Eq. (93), must also vanish at the sonic point if an accept- 
able solution is to be achieved. This does occur at N| = 0.375, and the cor- 
responding distributions of ^ , (j) , £ , and y are shown in the remaining 

four graphs of Fig. 15a. Figures 15b - e give similar results for Y = 2, 

2. 5, 4, and 6. To illustrate the effect of « , another calculation was done, 

at Y = 4, with (< =1.0, and cl - 1. The results, given in Fig. 16, differ 
from those found for |<^ = 10 chiefly by the fact that T is somewhat smaller, 

and Isl is closer to 2/5. This trend is qualitatively what would be expected; 

a larger value of /C strengthens the influence of the pressure gradient in 
drawing material laterally away from the axis of symmetry. The fact that 
more fluid is being extracted in this direction acts to retard the shock motion, 
i. e. , M is decreased. 

5. Comparison with the Symmetric Solution. The values of N/ found ir 
the above calculations were in all cases quite close to the value 2/5 that applies 
for the symmetric, constant-energy solution. Furthermore, the quantity y 
does not attain an appreciable value until some distance away from rj = 1 , 
where the density has fallen to a low value. Thus we might expect that, near j 



(96) 
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the shock at least, the asymmetric solutions will not differ greatly from the 
constant-energy solution. This is indeed the case. Figure 17, which compares 
the symmetric and centerline solutions for / = 3, shows that, along the center- 
line at least, the motion of most of the mass involved is well approximated by 
the solution for fsj =2/5. One may expect this trend to persist even for & 
greater than zero, suggesting that the Taylor solution will in general be an 
excellent approximation to the considerably-more-complicated asymmetric 
solution. The comparison shown in this figure is typical of the results found 
at other values of . This close similarity between the two solutions 

forms the justification for our use of the simpler constant-energy solution in 
making crater predictions. 

Furthermore, this close similarity appears to be a general feature. 
Admittedly, the evidence for this conclusion comes from a limited number of 
cases, in which a special pressure variation was used. Nonetheless, it is 
difficult to imagine how the true distribution of T(V}) could be radically dif- 
ferent from that used here. The function T (fj ) must always be zero at the 
shock, and must rise to a value comparable to the value of near the sonic 

point, i. e. , it must be of unit order. Future research in this area should 
examine other approximations for T , but it is considered highly unlikely 
that any contrary evidence will be found. 

So far as blast-wave theory is concerned, then, the energy of the pro- 
jectile plays the dominant role, its momentum being only of secondary impor- 

sjc 

tance. One plausible physical explanation is based on the experimental 

This explanation was suggested to the authors by Mr. Robert J. Yidal of this 
Laboratory 
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observation^^ that targets struck by hypervelocity projectiles often acquire 
momenta many times that of the projectile, implying that the material ejected 
from the target must also carry several times the projectile momentum. Thus 
it appears that the momentum of the projectile itself makes only a minor con- 
tribution to the over-all conservation process. 

A corollary of this conclusion is that the conditions of hypervelocity 

$ 

impact can be simulated by any experiment which duplicates the energy of 
the incident particle, irrespective of whether its momentum is correctly 
matched. In particular, any intense source of short-pulse electromagnetic 
radiation, such as the output of some currently available lasers, should be 
capable of providing such a simulation. Such an experimental technique 
appears to hold promise, and Appendix B discusses the basis for it in some 
detail. 

A final point must be made concerning the importance of energy versus 
momentum. One must not infer, from our use of a constant-energy solution, 
that "energy scaling" will hold, in the sense that crater volume is proportional 
to the projectile's kinetic energy over a range of impact speeds. The fact is 
that this theory predicts a result more nearly of the form R c V ^ 3 , due 
to the fact that y is allowed to vary with impact speed in an effort to fit 
the equation of state over the whole range. And equally important, the pen- 
etration law R c ^ v''\ sometimes described as "momentum scaling", does 

The fact that such a simulation is possible was first pointed out to the 
authors by Dr. Franklin K. Moore. 
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not imply that the projectile momentum is a controlling parameter. The 
theory presented here is derived from an analysis in which energy and 
momentum can both be conserved. Because the shock wave is such a strong 
agent in affecting the motion of the fluid particles, and because the low pres- 
sure region is, by comparison, so weak in this regard, the flow of the greater 
portion of the mass involved takes place as though it were caused by a 
spherically- symmetric disturbance. Given this fact, only the total energy 
need be considered. 

6. Isolation of the Sonic Point. This subsection describes the analytic 
means that are used to determine the exact location of the sonic point and to 
provide a valid solution at that point. In doing so, it is useful to work with 
the variables 


* = 


N <t> 


y = — 

0 rj* 




z = 


N -r 


rj o y 7 r) 

In terms of these variables, Eqs. (80), (87), (88), and (94) become (see 
Appendix A for details of the derivation) 


(97) 




(98) 


oU ~ -y 1) 

^ d 1 (*-N) Z -y 


dz 


Z M' 


£lf_ 

99 z 


(n>o) 


*-i e2 - z<2 - 


(100) 




- 3 - 2 - 


( 101 ) 


The vanishing of the denominator is now given by 



( 102 ) 


When this occurs, inspection of Eq. (99) shows that the solution will be single- 
valued only if 


at the same point. If (102) and (103) are both true at the same point, then 
(103) may be simplified to read 


H (fj ) , defined in Eq. (96), be zero. At the sonic point, where (102) and 
(103) are both true, the derivatives in Eqs. (98), (99), and (101) become 
indeterminate. The manipulations required to resolve this indeterminacy are 
given in Appendix C. There it is shown that the slopes at the sonic point 
(designated by an asterisk) are given by 




- 4 * (*"l) =0 


(104) 


which a little algebra will show to be equivalent to the condition that the functioi 




(105) 


where is a solution of the quadratic 


A *f £>s&i + C — 0 


( 106 ) 


and where 




N ) [z - (/- l) **-(3 r- l) <2 
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The functions A, B, and C appearing here are defined in Appendix C. They 
depend on the sonic -point values of , y. , z , -P , (f> , and 

T > in addition to their dependence on N and Y . Having found JL, and 
, the formulas presented in Appendix C permit one to find ctx./cCq'j* , 
f'T. <pT > T 1 ) , and -P ) , and with these, the solution can be continued 

smoothly through the sonic point. All that is required is the list of sonic - 
point values referred to above. These can be found by extrapolating from a 
set of solutions, at various values of h| , which approach more and more 
closely to the sonic point. The procedure is illustrated in Fig. 18; here are 
shown, in the y y plane, some of the numerical solutions obtained for the 
case Y =3, with X = 10, a, =1. The line marked "numerator = 0" is the 
locus of points at which Eq. (103) is satisfied, while that marked "denominator 
= 0" is the locus of points satisfying (102). By extrapolating these two loci, 
one may isolate the values and y* at which these two conditions are ful- 

filled simultaneously. Similar extrapolations permit a determination of 


1 )* > Y*. f 


, and r* • These, in turn, fix the values of all the slopes, 


so that the solution in the neighborhood of the singularity is complete. The 
pair of separatrices shown in Fig. 18 was found by just this method. 

Even with the solution determined in this fashion, there still remains 
some question as to the precise value of N to be identified with it. For 
practical application, there seems to be no point in specifying N/ to five- or 
six-figure accuracy. Thus, in carrying over results from the very fine scale 
of Fig. 18 to the "practical" scale of Fig. 15a, the whole solution is tagged 
with the number fsj = 0. 375. In actual fact, the solution shown in Fig. 1 con- 
sists of the solution common to<0.37500 and 0.37501 for /£> 0 . 055, plus the 
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extrapolation through the sonic point. 

7. One -Dimensional Solution. The remainder of this section presents 
a treatment of the disturbance field behind a strong planar blast wave. Becausi 
only one spatial dimension is present, a complete solution can be worked out. 
This solution illustrates a method for simultaneous energy and momentum con- 
servation in a self-similar solution. In addition, the formulas which are deriv< 
can be used to predict the craters that will be formed by the planar impact of 
slab-or disk-shaped projectiles. 

We make the perfect-gas approximation from the outset, choosing the 
value of 7 ' in each case according to the method outlined above. The problen 
is to determine the motion of a half- space of material resulting from the impul 
sive application of a large pressure to its free surface, which is then exposed 
to a perfect vacuum. The large impulsive pressure may be imagined to arise 
in various ways, for example by the detonation of an explosive at the free sur- 
face, or by the face-on impact of a disk-or slab-shaped projectile: 


EXPLOSIVE 

SLAB 




In the ensuing motion, a strong shock wave propagates into the target, while 


the free surface expands rapidly into the vacuum on the left 
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We may expect that a blast-wave solution will correctly describe the motion 
when the time after impact is large compared to the time during which the 
impulsive pressure was applied. Thus the free surface will lie, effectively, 
at # = — oo . The blast-wave solution of this problem has been discussed 
by a large number of authors. An excellent review of these is given 

by Mirels. 

The equations of motion are 

3 


3£- + 

o )t 


(f>U.) st 0 


(108) 


3 1 9 $, /° 



(109) 

(110) 


In a self- similar solution, the dependent variables are assumed to depend only 
on the similarity coordinate 



% 

*,(t) 


and are made dimensionless by the definitions 


** <t>(q) Pqi) 


i f(fj) 


(in) 


( 112 ) 


As before, all time dependence is eliminated from the differential equations if 
the shock advances as a power of the time 

* 5 =A t N (113) 

Thus the similarity equations become 

(<t>-q) f' * 4>'f (U4) 
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<P' - ±= ff~ <t> * y- = 0 (ns) 

- 2 f + (<P -i)< r ' - -^r (<P~i) v'=o <116) 


Solved for the derivatives, these become 

.±u± A L - *<*"») 


r = - 


N 

/-M 

N 


6 M) - 


yf 
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(0-)){(M) - it} 




si 


(Mf - f- 


with the usual boundary conditions at the shock 


<P0) = f (o = 


rw 


^(/) = 


/w 

/ 


(117) 


(118) 


(119) 


( 120 ) 


8. Determination of N . Before Eqs. (117) - (119) can be integrated, 

the parameter f\J must be specified. Just as in the axisymmetric case, the 

criterion that determines N is that the solution pass smoothly through the 

sonic point. The values of /sj that accomplish this feat have been reported 

in the literature for several values of / , and a few more were calculated in 

the present effort. These results are summarized in Table II, and are plotted 

in Fig. 19 for y between 1. 0 and 10. The distributions of pressure, density 

and velocity that are obtained in this type of solution are shown in Fig. 20 for 

(23) (27) 

the case / =7/5, which has the exact solution' ’ ' ' 

, - 5/2 c - 3/2 

f = <b(5 - 4iy) , $ = -^(2rf-l) , -P=^(5-4^) ( 121 ) 

In further analogy to the axisymmetric case, the values of N which 
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permit a smooth crossing of the sonic point are different from those cor- 
responding to conservation of the total energy or momentum (per unit area 
in the plane of the target surface). Denoting these respectively by g and 


, the appropriate integrals are 

5‘C*' (h = &***>[ (£r * -T t ty 2 ) * 'l < 122 > 

-o o 


n\ 

pu dt, - tytpdri (123) 

-oo J-O0 

Because ^ t , the energy integral will be independent of time only if 
N = 2/ 3, while the momentum integral requires N = 1/2. Neither of these 
values, however, leads to a physically acceptable solution. Instead, the 
solution for N = 2/3 is symmetric about oL = 0, while that for N = 1/2 becomes 


double- valued. The solution passing through the sonic point is the only one 
that starts at the shock, and yields a velocity distribution like that of Fig. 20, 
in which material moves rapidly to the left at large negative values of 

The values of N corresponding to constant energy and constant 


momentum are compared with the curve of N(y ) in Fig. 19. It should be 
noted that for / in the range from 2.0 to 10, N(/) is relatively constant, at 
a value not much different from the value that holds for constant energy. As 
a corollary of this, we find once again that the constant- energy distributions 
of pressure, density and velocity are very nearly the same as those having 
the smooth sonic-line crossing. Figure 20 illustrates this point for the 


case =1.4. (The constant-energy solution used here is described in 

detail below. ) The close identity between these two types of solution is 
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exactly the same behavior as found above in our approximate treatment of 
the centerline equations. The fact that the one-dimensional results are 
found without approximation (there is no term analogous to Y ) may be 
taken as evidence for the generality of the conclusions reached in the axisym- 
metric case. 

9. Simultaneous Conservation of Energy and Momentum. There can 
hardly be any question about the fact that the solutions which cross the sonic 
line smoothly are the only acceptable ones. Having chosen N so as to 
achieve such a smooth crossing, however, it is now impossible to satisfy 
either of the conservation conditions, in the form given by Eqs. (122) and 

25 

(123). One means of resolving this difficulty is pointed out by Zeldovich. 

He suggests that a small portion of the mass be ignored in calculating the 
total energy and momentum. This small mass, which originally lies at the 
free surface and is strongly compressed during the initial stages of the 
impact, acquires an entropy during the impact phase which is different from 
that predicted by the similarity solution. This entropy depends only on con- 
ditions at the shock, through integration of the differential equations, and 
has nothing to do with the details of how the shock was formed. The entropy 
of each particle remains constant, once it has been processed by the shock. 
Thus the small portion of mass processed during the impact phase always 
bears the imprint, so to speak, of this phase, and its motion is never cor- 
rectly given by the similarity solution. Zeldovich’ s argument is that in 
seeking a self-similar solution, one should consider only that mass whose 
motion is expected to be properly described by such a solution. His argument 
derives considerable support from the fact that non- similar numerical 
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( 27 ) 

solutions' of the partial differential equations in ^ and £ fair 
smoothly into his similarity solution when the time is large compared to 
the duration of the impulsive pressure. 

Mathematically speaking, the exclusion of this mass corresponds to 
the replacement of the infinite lower limit in Eqs. (122) and (123) by a 
finite number — (£) . The introduction of a second free parameter per- 

mits two conservation conditions to be satisfied. In the remainder of this 
section, we work out in detail the relations that follow from neglecting this 
small amount of mass. 

The point - 4 0 (t) ( 4 ^ > 0) , which forms the boundary of the neglected 

mass, will be far to the left of the origin, since the mass neglected is small. 
As a prelude to modifying the conservation integrals, then, we must first 
determine the asymptotic form of the solutionat large ne gative' value s of Y| 

The key to finding such an asymptotic solution is to neglect the pressure- 
gradient terms in the momentum equation, Eq. (115). We have seen before 
that in these solutions the particle velocity is large at large negative r| , 
with a slope of order one. The pressure, on the other hand, is approaching 
zero, signifying that its derivative is also becoming extremely small. Neg- 
lecting this term yields the simple solution 




N 


(124) 


With (j> known, the density ^ can now be found from Eq. (114) as 

t*A, hf 7 ^ ,125) 

Finally, use of <f> and ip in (116) gives the asymptotic pressure formula 


a- 


t * A. 


I-N 


(126) 


RM-1 655-M-4 


57 


Substituting back into Eqs. (11 4)-( 1 1 6) shows that these approximations 
are internally consistent. 

( 29 ) 

The constants and are not independent. Lees and Kubota, 
among others, have shown that the similarity equations in general have the 
solution 


f-Z 


l-N 

N 





(127) 


where, from conditions at ^ =1, K, is found to be 

/ N 7>/ 

K - (i±j! 

N l “ n A. ,\r 


zO'-D 


(128) 


Equation (127) is derived by dividing (114) by < an{ l sub- 

tracting from this the result of dividing (11 6) by •£ (<J b~tj) • The relationship 
between Aj and is found by substituting the asymptotic pressure and density 
solutions into the general integral, Eq. (127). One finds that 


A* = 


20-M) 




(129) 


The quantity A^ is found by integrating the similarity equations out to a large 
negative value of fj , where the asymptotic behavior is reached. For 
example, in the case / ® 1.4, the exact solution given in Eq. (121) shows 
that A, = 3/16. In general, A, depends on / . 

If we now neglect a small amount of mass (mass per unit area), we 


have 




r '*o& pi* 

= ( y JL = fi. At” tdri 

J-oc 


(130) 
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The integral here can be evaluated directly. By writing the continuity 
equation in the form 


and integrating by parts, we find in general that 

ip drj = - (Q - 7 )] 

U CL/ 

Thus 

£ **1 = = -ir 


( 131 ) 


(132) 


(133) 


and 


At" 


/-N 

M 



(134) 


In order to make 77} independent of time, we must choose 




I— N 


k 2 > 0 


(135) 


which gives 


N 


».=A*,K i -4/*- ft, 


(136) 


The parameter is the second free parameter which, with A, will permit 
two conservation conditions to be satisifed. 

If we neglect the mass 7T( 0 in the energy-conservation integral, we have 


® = Po *s V 


r 1 




(137) 
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This integral, evaluated at the upper limit, yields a constant. At the 
lower limit, we may use the asymptotic solution. Doing so reveals that 
the internal energy contribution is negligible compared with the kinetic 
energy term, and we have 


A N t 


y^llf "" df l ■+ const ♦o(|nj* 

y -Vi, J 


(138) 


The second term in the curly brackets corjftes from the upper limit, the third 
term from the internal energy. Carrying out the integration gives 


*-pA*** 1 *. h 


.3-r 


I— N 


3 - 


+ const. + o(|rj ( 


3-TT 


l-N 


<-N 


= & 2 * 


A 3 -(2-3N) f ,_ M 


| A ' 2-3N 


K * ,_N t ' W +■ const. + 0 (t 3 i 3W )| 


(139) 


2-3N 

^ o A 3 A ' IzN ^ “W- 




2 2-3N ’'a. 

as t becomes large (this is the limit at which we expect the blast-wave 
solution to hold), the second and third terms here become negligible, since 
N is less than 2/3, and / greater than one. 

Neglecting the same amount of mass in the momentum integral gives 

(140) 

This integral can be evaluated exactly. Re-writing Eq. (115) as 

zh 

N 




--TT 


(141) 
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and integrating, we have 

rb 


[ ~ N 

N 


V<pdrj-hfj + J” f (<p-rj) drj = 0 


Consider only the last term. We've already seen in Eq. (132) that 
An integration by parts gives 

? (4>- r j) <t>' dr/ = W- 7 )^] + f ^ 

Jd, * 

Substituting this in (142) gives the desired result 

J 1o 


(142) 


(143) 


(144) 


(145) 


The right-hand side is found to be zero at the upper limit, and at the lower 
limit, the pressure contribution is small compared with the second term. 
Thus 


I yfdtj x 
\ 



I — N 
2N-I 



(146) 


As |f | 0 |— >00 , the integral vanishes. Thus, if we were to include all the mass, 
and describe its motion by the self- similar solution, we would find zero net 
momentum. In fact, Zeldovich takes this point of view, and apparently con- 
siders the solution valid only for such a momentum condition. However, it 
would appear to be more consistent to neglect mass in the momentum integral 
as well. Doing so, we have 




r 


= p 0 At* NAt*"' 


Ai(l-N) ^ ~i fi"~ , 1-2N 
H (ZN-0 ** Z 


= e 0 a 2 a, 


l-N 

2N-I 


I- 2 N 


X 


l-N 

2 


( 147 ) 
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This relation, together with (139), enables us to find A and in terms of % 
and p . Note that the velocity of the impacting slab is related to these 
parameters by 


zt __ 2N-I 
p “ Z- 3N AK * 


(148) 


while the neglected mass is of the same order as that of the projectile 


? 2 _ „ A . 0-N)(2-3N) 

(ZN-i) x 


-P77 N(2- BN/) 

“ (ZN-IP 7770 


(149) 


Solving for A and K a in terms of and P gives 


K 2 - 


_ ) 4-A 





\ l-N 


J 


A - G(?) 




where 


&(/;» 


ZN-I 


( 2 / S /-/) 3 

4A . 

\ 1 


\ l-N 
r 

J 


(150) 


(151) 


(152) 


The actual solutions of the similarity equations have never been reported, 
except for ^ =1.4. Thus the variation of A with is not presently known. 
When this information becomes available, it will be possible to compare the 
final formula for the shock propagation 
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( 153 ) 


l s (t) = Gr(t) 


? 


2 




with the one -dimensional constant- energy result 

3.V/ 3 


' Ur jit) 


Po 


(154) 


Figure 21 shows the function 1 it) > found by a hand computation using the 
exact solution which is given, for example, in Ref. 6. The constant- energy 
distributions of pressure, velocity, and density are compared with the 
asymmetric result in Fig. 20 for Y -1.4, 

The "craters” corresponding to these shock-wave time -hi stories may 
be found as the shock depth at the instant when the pressure behind the shock 
has decayed to the level IP . Denoting this quantity by , a little 

algebra gives, for the asymmetric solution 


i 


c 


Git) NV 

Po £ \ Y+t 


a 



- > 
✓ 


(155) 


and for the one-dimensional, constant-energy solution 




2 J_ 

(7+i) J- 0 (t) P 


(156) 


The method described above can be applied, with analogous results, 
to the axisymmetric case. Far from the impact point, we may imagine that 
the flow is confined within a cone for which Q & Tf 
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Within this cone, all quantities have a negligible variation with 0 , and 

the velocity component cO is small compared with ^ . The problem is 

then essentially the same as that treated above, and the same series of steps 
can be taken. The only new term present in the axisymmetric case is the term 
in the continuity equation, which merely has the effect of changing cer- 
tain exponents. The conclusion remains the same, namely that simultaneous 
conservation of energy and momentum can be achieved by neglecting a small 
amount of mass. 
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II. BOUNDARIES FOR APPLICATION OF 
THE BLAST WAVE ANALYSIS 


A . The Space Vehicle Environment 

The probability of collisions between space vehicles and meteoroids 
of sufficient size to cause serious damage is too small and too variable to 
allow convenient verification of empirical relationships by satellite experi- 
ments. Therefore, it is necessary to improve prediction methods as much 
as possible. Diverse sources of evidence can be used. Among these 
sources are theoretical analyses, hypervelocity impact experiments, studies 
of meteors, meteorites, meteor craters, and impacts of smaller meteoroids 
on satellites, etc. A survey of information pertinent to this problem has 
been published by Davison and Winslow. 

There is little agreement concerning the properties of the impacting 
particles. Density estimates and assumptions range from 0.05 gm/cm^ 
to 7.8 gm/cmV^ Dense meteoroids are assumed to be of asteroidal origin 
and low density meteoroids are assumed to come from comets. Whipple^^ 
indicates that in space about 10% of sporadic meteoroids and not more than 
20% of all meteoroids are of asteroidal origin. Of the meteoroids observed 
in the vicinity of the earth about 50% are sporadic and 50% are in streams. 
Relative velocities range up to 80 kilometers per second (assuming a par- 
ticle in a closed orbit about the sun moving at 42 kilometers per second 
approaching the earth in a direction opposite to the earth's motion at 30 
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is 28 km/sec. The average velocity of smaller meteoroids is less. 9 

Whipple* 32 ^ chooses a velocity of 1 5 km/sec for the smallest meteoroids 

and assumes an arbitrary variation of velocity with magnitude. The fainter 

meteoroids and those detected by radar observations tend to be more spora- 

(35) 

die than the brighter meteoroids. ' The mass of individual meteoroids 
has been determined from brightness measurements on the Harvard Photo- 
graphic Meteor Program. Meteors of visual magnitude zero were found to 
have masses of the order of 25 grams. Based upon previous rough esti- 
mates^ 32 ^ of the penetration of meteoroids through structural materials, 

-2 -4 

the mass range of interest is from about 10 to 1 0 grams. Conversion 
of visual magnitude to estimated mass is accomplished by the following 
formula* 34) ^ 

M ~ 2.5 Anfrio 

in which M is the magnitude 

1Y1 0 is the mass of a zero magnitude meteor 
and 7)f is the mass of the meteor 

Photographic techniques are used for meteoroids down to the 5th magnitude. 

(32) 

The use of radio echo techniques may be extended to the 12th magnitude. 

A comparison with these mass estimates indicates that the most direct 
evidence of interest for this penetration problem will come from radio 
echo measurements. 

(37) 

In a study of photographic meteor wakes and trains McCrosky 

found that the particles which break off to cause the wake, have masses 
-5 -6 

of 10 or 10 grams. This evidence of the mode of fragmentation and the 
size of fragments making up the larger meteoroids lends additional weight 
to Whipple* s suggestion that most of the meteoroids are low density bodies. 
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The original estimates of the density of photographic meteoroids have been 

(38) 

revised upward recently. On the basis of new estimates of luminous 

efficiency a one gram meteoroid having a velocity of 30 km /sec. has a density 
of 0.44 gm/cm . These estimates result from experiments in which par- 
ticles of known properties were projected down through the atmosphere and 
simultaneously photographed from the earth. 

B . Energy Transformation During Impact 

An attempt has been made to present the hypervelocity impact pro- 
blem as a chronological sequence of events involving realistic solids. This 
approach leads to consequences that are difficult to unify in mathematical 
terms but which, it is hoped, will give some insight into the nature of im- 
pact processes. 

During the brief period of meteoroid impact there are several stages 
of energy transformation. It is important to consider the time history of 
these stages because target characteristics will partially determine which 
stage predominates in the production of damage to the space vehicle. Both 
the meteoroid and the target may undergo polymorphic phase changes, how- 
ever, the time occurrence of such phase transitions will not be simultaneous 
unless the materials involved are identical. It is extremely unlikely that the 
meteoroid will be identical in mechanical properties with the satellite target 
materials. 

Table 3 lists the stages of energy transformation undergone by the 
meteoroid-target interface under hypervelocity impact conditions. Every 
impact will not necessarily involve each of the listed steps since the impact 
kinetic energy varies over a very wide range. 
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The temperatures attained in satellite materials during impact are 
important for several reasons. The temperature attained is a measure of 


the energy absorbed into thermal motion. It also affects the assumption 
of fluid properties in the blast wave analysis and the ejection of material 
from the cavity in the final stages of the impact process during which one 
would expect that fluids would be lost from the cavity but solids would be 
retained. 

Solids may be heated or cooled by impact. The dominant factor 
affecting the sign of the temperature change is the Gruneisen factor (Pj) 
The Gruneisen factor for the ^-th mode of a particular lattice vibration is 

n 

f" V 

so that f\- is a measure of the change in lattice vibration frequency for 

t 

a given change in volume. The role of the Gruneisen factor in the equation 

(39) 

of state of solids is indicated by Slater' ' in a derivation based upon the 
statistical mechanics of a system of oscillators. The equation for the 
pressure is 


A= A, r i 




in which the first term on the right is the pressure at absolute zero and the 
summation is the so-called thermal pressure. 7? is Planck's constant 
and lc is the Boltzmann constant. 

The role of the Gruneisen factor has also been indicated by 
Eastabrook^^) starting with an equation derived using the Euler reciprocity 
relation 

_ (2£_j 

3T/ V VsvAr 
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in which T is the absolute temperature and S is the entropy. This may be 

VrN 


rewritten using Maxwell's relations 


so that 


-&U U-M 


3V, 


p — V (C T 

c v/ 


(157) 


in which o£ v is the volume coefficient of thermal expansion = 
is the isothermal bulk modulus = — V 
and C is the heat capacity at constant volume = T Hr~j 

P^ Gruneisen ’’constant 1 ' which is an average of the 

taken in a way so that the equation is satisfied. Since the various P . > may be 

f 5 

either positive or negative, wide variations in the values of P necessary to 
satisfy Eq. (157) are observed. The factors C y , V and are positive 
quantities, so that P will be negative when the volume thermal expansion co- 
efficient is negative. Some materials have negative values of thermal expansion 
coefficient, at least over part of the temperature range for which these measure- 
ments have been made. Among these materials are vitreous silica (SiO^), silicon, 
germanium, diamond (C), indium antimonide (InSb), several lithium aluminum 
silicates and uranium pyrophosphate (UP^O^). The thermal expansion coefficient 
of uranium pyrophosphate is of special interest because it is negative at high 
temperatures . 

Temperatures during impact have been calculated by several investiga- 
tors, Wide variations in temperature are found depending upon material 
properties. The calculations by Walsh et al^ ^ of temperature rise along an 
adiabat are obtained by integrating the equation 
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r T £a r v rM 

\ T ""l v 


V 


This calculation depends upon knowledge of the process and an estimate of 

the volume dependence of P . 

(41 ) 

Benedek's' ' temperature calculations depend upon knowledge of 

(42) 

the volume dependence of the cohesive energy. The equation used is 

V(/ 1 ) -v Co) = E (yu) - B(0) - $ (/0 - § (0) 


in which 


_A c 


2 2 


and y6( = volume compression 

V = vibrational energy 
E = total internal energy 
^ = cohesive energy 

and C and S are the constants in the equation 

Ur - C + Sa, for the shock and particle velocities 

Knowing the vibrational energy, the temperature is calculated by use of the 

Debye approximation for the specific heat. 

(43) 

Wackerle 1 s' ! calculations are based upon the following formula: 


T" — T** « 


<k) 


n»/* r>^ 

e e 




in which is f 1 evaluated from the STP values of od v » Cp and 

adiabatic bulk sound speed, a,nd and T H are the initial and final temper- 

atures along the Hugoniot. The results of the calculations described above 
indicate roughly the expected dependence upon thermal expansion coefficient 
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and P ; that is, the low expansion materials show, initially at least, 
smaller temperature changes than high expansion materials. The calculated 
temperature rise of vitreous silica is only 5°C at impact pressures up to 
262 kilobars. 

It is evident that the space vehicle designer can have a degree of 
control over the temperatures attained during impact. In regions in which 
it is desirable to absorb the maximum amount of energy in thermal vibra- 
tions, the choice of materials with high P is indicated. In regions in 
which it is desirable to prevent melting in order to keep material in the 
crater after impact one would choose materials with low values of P 
Therefore, one expected feature of these configurations is a sandwich struc- 
ture with the high P material on the outside and low P material on the 
inside. 


C . Stress Wave Velocities in Solids 

1 . Elastic Stress Waves . At low velocities the impact of a small particle 
upon the outer (extended) surface of a satellite will introduce a stress wave 


in the target which involves the dilatational wave velocity given by 

r . * \ n'A 


- 

_< 


Y (hjO 

fd + JXi-zs) 


(158) 


where Y = Young's modulus 

i) = Poisson's ratio 

Note that the bulk compressional wave velocity (Eq. (158) ) is significantly 

(441 

different from the compressional wave velocity found in a long thin bar' ' 


RM-1 655-M-4 


71 



(i.e. , where the bar cross section is small compared with a wave length). 

The change in velocity which occurs as the wavelength approaches the cross- 
sectional dimension is due to the increase in lateral restraints. For the 
high frequencies associated with short duration meteoroid impact, the higher 
dilatational wave velocity is more appropriate. In isotropic solids, Poisson's 
ratio is always less than one-half and greater than zero, hence Q ^ C c 
(the velocity in a long thin bar). 

For isotropic solids, the fastest elastic wave travels with the bulk 
dilatational wave velocity. Hence, meteoroid impact will be in the hyper- 
velocity region if the impact velocity is greater than the fastest elastic wave 
velocity in the target material given by Eq. (158). 

In Table 4, the bulk dilatational wave velocity is tabulated for several 
plastics, metals, ceramics and glasses. Beryllium has the highest stress 
wave velocity among the metallic elements. Corundum ( oc - alumina) and 
silicon carbide are among the highest for ceramic bodies. Therefore, im- 
pact velocities above 13 kilometers per second are probably in the hyper- 
velocity region for elastic wave propagation in all materials. 

The transition region between subsonic and supersonic velocities is 
defined here as the velocity range lying above the lowest elastic wave 
velocity (Rayleigh waves for short duration impact) and below the highest 
elastic wave velocity (bulk dilatational wave velocity). The transition 
region will become broader as Poisson's ratio increases towards 0.5 so 
that different target materials will have both different critical impact velo- 
cities (i.e. , where V— C, ) and different transition velocity ranges. 

In compression, most materials exhibit an increasing Young's modulus 
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with increasing strain, i.e. , 


d'G't y {_ d'Cr c _ 


(159) 


ds 1 d.& 

Similarly, the bulk modulus ( ) usually increases under isostatic compres- 

sion. In isotropic materials, the relationship between the shear modulus 
( G- ), Young's modulus and bulk modulus is given by 

AKY 


& = 


9l(Y 


( 160 ) 


Since both y and % increase under isostatic compression, the shear 
modulus of isotropic solids will also increase under the same conditions. 

In some solids the elastic moduli increase with temperature. In these cases 
isostatic compression will cause a decrease in the elastic moduli. Negative 
pressure coefficients of elastic constants have been measured at room temper- 
ature for fused silica and Pyrex glass both of which have positive temperature 
coefficients of the elastic constants below 500° C. 

The behavior of the elastic constants of single crystals under isostatic 

(45) 

pressure is more complex. The shear stiffness coefficient ( ) de- 

creases with pressure in certain cubic crystals, notably KC1 and RbCl.^^ 

In the less symmetric crystals, several of the elastic constants may have 
negative pressure coefficients. Alpha-quartz, for example, has at least two 
negative pressure coefficients as shown in Table 5. In alpha-quartz, the 
shear stiffness along the optic axis ( Q ) decreases with pressure while 
the shear stiffness along the x and y- axes ( C ^ ) increases with pressure. 
Clearly, generalizations about the effects of pressure upon the shear modu- 
lus based upon atomic force-constant considerations must include more than 
the nearest neighbor interactions if negative pressure coefficients are to be 
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predicted. 

2. Plastic Stress Waves . The elastic wave velocity relations given 

in the previous section are valid for small amplitude stress waves. With 

increasing stress-wave intensities, the velocity relations become nonlinear 

functions. These nonlinearities arise from two different factors. The first 

is the nonlinearity caused by higher order terms in the stress- strain relations 

(471 

for finite deformations. The theory for finite amplitude stress waves is 

a complete subject in itself and will not be treated here. The second is the 
nonlinearity associated with plastic deformation where the tangent modulus 

associated with infinitesimal amplitude stress 
waves is dependent upon the strain. High-intensity stress waves whose 
stress amplitude exceeds the elastic limit are known as plastic waves. In 
general, the plastic -wave velocity is dependent upon both the strain and the 
strain rate.* 48 * 

Strain-rate effects in metals have received extensive experimental 

study. Each class of metals or alloys appears to behave in a different 

(49) 

manner. Dislocation theory' 1 provides a theoretical explanation of strain- 
rate effects. The evidence suggests that for crystalline materials, strain 
rate dependence is associated with motions of dislocations. Dislocations 
are arrays of lattice-point defects which may take several different forms, 
e.g. , a line, a loop, a helix, a screw or combinations thereof. Each type 
of dislocation has a characteristic velocity which is dependent upon the host 
lattice structure and the absolute temperature. Hence a general discussion 
of strain-rate dependence becomes quite complex. The possibility of dis- 
location motion at velocities above the highest elastic stress wave velocity 


of elasticity 
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is the subject of considerable controversy. If the velocity of any disloca- 
tion motion may not exceed the fastest stress wave velocity ( ), under 

hypervelocity impact conditions, dislocations ahead of the penetrating pro- 
jectile are not free to move and the intrinsic strength or dynamic yield 
point is determined by the interatomic forces alone. 

Plastic waves have been studied extensively under tension impact 
conditions. Because of the decreasing slope of with increasing 

strain for most metals, the plastic wave portion of a tension impact distur- 
bance travels with a lower velocity than the elastic wave portion causing the 
stress pulse to become elongated at positions remote from the point of dis- 
turbance. 

3. Shock Waves in Solids . Since the elastic moduli for isotropic solids 
usually increase under compressive stress, the several stress wave velocities 
will also increase with compressive stress. Thus, attempts to propagate in- 
tense stress waves will result in the generation of a discontinuous wavefront. 
Stress waves which require discontinuous functions to describe the wavefront 
are known as shock waves. The shock wave velocity will usually be greater 
than the corresponding elastic wave velocity in a given isotropic solid. 

In a gaseous medium, the thickness ( £ ) of an intense shock front is 

equal in length to a few mean free paths of the gas molecule s . By 

analogy, the shock front thickness in a solid medium is on the order of a 
few phonon mean free paths. In gases, thermodynamic equilibrium is reached 

❖ 

The term shock wave as used in the literature generally refers to an intense 
bulk-dilatational wave. However, in the broader context used here, shear- 
mode shockwaves are also possible. 
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only at finite distances behind the shock front. Similarly, thermodynamic 
equilibrium in solids behind an intense shock front requires distances 
equivalent to several phonon mean free paths. The lack of thermodynamic 
equilibrium within the shock front can be illustrated by comparison of the 
impact velocity with the instantaneous velocity of lattice vibrations (near room 
temperature). Assuming, for convenience, that the lattice vihrations behave 
as harmonic oscillators, (i.e. , ) then the instantaneous lattice- 

point velocity will be given by 

= cJ X oO t 

where Ct) = Zir-f 


f = frequency of lattice vibrations 

X = maximum displacement of the vibrating atom within the 
lattice unit cell 

For a maximum displacement of one-tenth thelattice spacing ( ^ ), a lattice 

13 

frequency of 10 cycles per second (infrared band), and a nominal lattice 
spacing of 5 angstroms, the maximum instantaneous lattice -point velocity is 



Thus, the maximum lattice-point velocity is less than the impact velocity 
in the hypervelocity region. This simple calculation is sufficient to illustrate 
that even on an atomic scale, impact velocities in the supersonic range repre- 
sent a gross disturbance to the lattice point displacements at the impact 
interface. 

- 8 

The phonon mean free path ( A ) in alpha-quartz at 0°C is 40 x 10 cm 
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(51) 

or about 8 lattice spacings' ' given by 

A - 

where = thermal conductivity 

H = heat capacity per unit volume 

C = "mean sound velocity” - space averaged elastic wave velocity 
The shock front thickness ( f 5 ) in alpha-quartz would be in the order of one 
micron near room temperature ( t >A). 


D. Collisions Involving Porous Bodies 

The calculations involving porous meteoroids are especially important 
because of the relatively great frequency of occurrence of cometary meteoroids. 
Impacts and shock waves in porous materials have been studied to a very 
limited extent. Altshuler et al obtained shockwave data for porous 

(specific gravity = 5.52) and non-porous (specific gravity = 7.85) iron to 
5 x 10^ atm. The equation of state for tuff, a porous rock, is also available . ^ 

In the following sections, the conditions during impacts between porous 
materials and non-porous materials are compared to indicate the differences 
caused by the presence of pores. 


Case I - Impacts of Porous Meteoroids Against Solid Surfaces 

Calculations of the shock-wave characteristics for impacts between 
porous and solid iron, using Altshuler's data and assuming planar shock waves, 
will be used to indicate the effect of porosity. The Hugoniots are shown in 
Figure 25. The presence of porosity results in an increase in the pressure 
and energy per unit mass at a given density of material in the shock front. 

Since the pressures in the meteoroid and space vehicle surfaces must be 
equal at the point of impact, the shock and particle velocities can be deter- 
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mined. For example, consider the case of a porous iron meteoroid colliding 
with a solid iron space vehicle surface at an impact velocity of 9 km/sec. 

If we ask which combination of particle velocities for the porous and non- 
porous materials add up to the impact velocity, we find using Figure 26 
that at equal interface pressure the particle velocity ( LL, ) in the solid vehicle 
skin is 3.8 km/sec. While the particle velocity ( OL, ) in the porous "meteoroid" 
is 5.2 km/sec. The shock velocities and density ratios can be read from 
Figure 27. The energy per unit mass can be calculated from 



The results are given in Table 6. The energy per unit mass in the porous 
"meteoroid" (1.4 x 1(3^ ergs/gm).is significantly greater than that in the solid 
vehicle skin (0.76 x 10^ ergs/gm). It seems likely that this energy distri- 
bution will result in decreased meteoroid penetration due to greater evaporation 
which in turn is due to greater conversion of kinetic energy into thermal energy 
at impact. This conversion into thermal energy is illustrated in Figure 28, 
which has been revised from Altshuler's paper. Curve /jr-^ is the schema- 
tic shock compression curve for a non-porous material. The thermal energy 
( ) added to a unit mass compressed to /lAJ is indicated by the area of 

triangle ft f ft less the area of triangle which is work required to 

compress the material to /it at 0°K ( £.^ ). The same type of material 

(i.e. , iron) containing pores will be compressed by shock waves to various 
states along nr^ft^ If we neglect the small amount of energy required to re- 
move the pores by compression at absolute zero, (indicated by the dotted 
line), the extra thermal energy added (AE T ) to the porous material to com- 
press it to AT, is given by the area ft ft z ' >r to ' v ofti • Presumably, the 
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increased temperature and pressure at will cause increased evaporation 

or melting of the porous material when the pressure is released. 

Cases II-V 

When a non-porous meteoroid collides with a porous vehicle skin 
(Case II), the particle velocity in the meteoroid is less than that in the vehicle 
skin (see Table 5). The shock velocities and density ratios are also reversed 
and the sign of the internal energy difference changes. When a non-porous 
meteoroid collides with a non-porous vehicle skin (Case III) at the same velocity 
as the preceding cases, the resulting interface pressure is increased and the 
internal energy increase per unit mass is the same in both bodies. When a 
porous "meteoroid" collides with a porous vehicle skin (Case IV), the initial 
interface pressure is lower than in the previous cases and the internal energy 
increase per unit mass is the same in both bodies and is roughly the same as 
the value in Case III. To increase the impact velocity so that the initial inter- 
face pressure is the same as that in Case III with porous materials requires 
an impact velocity of 12.2 Km/sec. In this case (Case V) the particle and 
shock velocities and the density ratios are increased. The internal energy 
increases per unit mass in this case is 2.0 x 10^ ergs/gm. 

It is evident that even a small amount of porosity significantly effects 
the shock variables. The most important of these effects is the increased 
thermal energy per unit mass. In the initial phase of crater formation, the 
increased energy converted to thermal motion becomes unavailable to do 
structural damage. 

If the materials are very porous, as seems to be the case for cometary 
meteoroids, the impact area is increased considerably and the initial interface 
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pressure will be decreased. Due to the large increase in impact area it may 
be better to calculate the crater depth using a one -dimensional theory. 

E. Crater Wall Criterion 

1. The Intrinsic Strength of Crystals. At low strain rates the shear 
strength in crystalline materials is determined by the type and distribution 
of dislocations. At high strain rates or when the geometry hinders dislocation 
motion, the measured strength approaches the theoretical limit imposed by 
interatomic forces. 

Under isostatic compression, the engineering strength of most non- 

(52) 

porous materials increases markedly. Bridgeman's data on static 

compressive strength of several metals and ceramics for pressures in the 
region of 400, 000 psi are given in Table 7. The static compressive strength 
increases with hydrostatic pressure becoming manyfold higher in normally 
hard and brittle materials. These data provide evidence of the importance of 
relaxation processes to static and dynamic mechanical properties. The pro- 
cesses are highly restrained under isostatic pressure so that the static 
strength approaches the dynamic strength. 

A method, originally due to Frenkel, for estimating the theoretical shear 

(51 ) 

strength of crystals is described by Kittel. If one considers the stress due 

to displacement of one perfect plane of atoms relative to another as shown in 
the sketch below, the stress to a first approximation can be represented by a 
sine function 



in which O' is the stress, G- is the shear modulus, a, is the space between 
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atoms in a plane, d is the spacing between planes and X is the displacement 



As atom moves over atom (jy) the stress increases, passes through a 

maximum and then decreases to zero. When atom © is directly above atom 
® » it is in a state of unstable equilibrium with respect to forces in the X 
direction. If the stress at the maximum value of this function is considered 
to be the critical shear stress, the value of the critical shear stress ( Qr^ ) is 
given by 0£. = &< 2 , ^ 2 It'd so that when a, and d are approximately equal, 

as for the case of shear in a K, 100} direction on a {/00} plane in a cubic 
crystal 

^ -fir 


in which G- is the shear modulus given . This approach to calculation 
of the critical shear strength has been refined, according to Kittel/ 51 ^ by 
J. D. Mackenzie using an estimate of the form of the interatomic forces and 
by "consideration of other configurations of mechanical stability that the 
lattice may develop as it is sheared. " These effects reduce the critical shear 
strength to about G/30 in some cases. 

Strengths as high as those predicted by this atomic model have been 
observed for a number of materials in whisker form. In these crystals the 
geometry restricts dislocation motion. Hoffman^ ^ lists the Young's modulus, 
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ultimate tensile stress, and the tensile stress to Young's modulus ratio for 
the strongest whiskers prepared from a number of materials. Several mater- 
ials from this list are retabulated below with addition of the dynamic shear 

modulus and the maximum shear stress (o' ) to shear modulus ratio 

s mux 

(assuming 01 = - 5 - O’. ). 

s max c 

COMPARISON OF THEORETICAL AND OBSERVED STRENGTH 



Y oung 1 s 
Modulus 

Dynamic 

Shear 

Modulus 

Ultimate 

Tensile 

Stress 


max 


(E) 

(G) 

( 0 * ) 

E 

& 

Iron 

29 x 1 0^ psi 

11 .9 x 10 6 psi 

1 . 9 X 10 6 

1/15 

1 /1 2 

Copper 

18 

6. 7 

0.43 

1 /42 

1 /32 

Sapphire 

(corundum) 

74 

23. 5 

1 . 7 

1 /43 

1 / 28 

Carbon 

1 x 10 6 


0.88 

1/11 

— 


The ratios <T s /q, can be compared with 1/2 V and 1/30. The observed 
ratios correspond quite well with the theoretically derived values! Since the 
observed strength of whiskers increases with decreasing whisker diameter 
and since the whiskers may have failed in tension rather than shear, the re- 
sults for smaller diameter whiskers failing in shear are expected to yield 
larger ratios. The possibility of still larger strengths is confirmed by the 
recent results of McQueen and Marsh v who showed that the yield strength 
of copper is in excess of 150 kilobars for shockwave experiments." If the 
maximum shear stress is one half the tensile stress at yielding, this indicates 
a shear stress of 1 . 02 x 1 0^ psi. Therefore, for copper & s /G- = 0.152 which 
can be compared with 1/2 Tf = 0.159, giving an unexpectedly close comparison. 
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McQueen and Marsh also indicate a critical pressure on the Hugoniot which 
corresponds to the yield point {oCP/cLV - 0). The negative pressure at this 
point is a theoretical tensile strength. These theoretical strengths, with 
corresponding shear stresses and other data are given in Table 8. The cal- 
culated values of based on this information are somewhat larger than 

1/2 Tf . Based upon all of the information given above, the critical shear 
strength = G/2f Y is considered to be the best available value for the shear 
strength at high strain rates. 

2. The Influence of Shear Strength on Crater Formation . It has been a 
common practice in fitting impact equations to experimental data to normalize 
the impact velocities to the velocity of sound in the undisturbed medium. The 
presence of a transition region, the velocity range near the sound velocity in 
which there is a change in the relationship between crater depth and impact 
velocity, has served as a partial justification for this practice. However, the 
reason for the change in the relationship between crater depth and impact 
velocity has not been clear. It seems reasonable that the processes governing 
failure (i. e. , crater formation) should depend primarily on the properties of 
the disturbed medium rather than the undisturbed medium since there is a sub- 
stantial pressure discontinuity at the shock front even after the crater has 
reached its final form. 

Pressure profiles for impact of iron on iron and a natural meteoroid on 
tuff have been calculated numerically by Bjork.^’^ These profiles indicate, 
as might be expected, that the pressure gradients in the disturbed material 
are greatest near the crater lip. It seems likely that shearing of the material 
at the crater lip is the first requirement for ejection of material from the crater. 
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at least in regions far from the impact point. 

If one considers a cube of material at the crater lip, and assumes that the 
material on the inside (toward the impact point) has already been removed or 
melted, it is possible to visualize the shear stresses on this material. The 
cube is pressed from behind by the high pressure material in the cavity. Motion 
is resisted by the pressure of the material moving ahead of it and by inertial 
forces and shear stresses. As the crater wall reaches its final position, the 
motion of material slows down and stops while the pressure gradients and shear 
stresses remain substantial so that, at this point, the inertial forces in the 
tangential direction can be neglected. Therefore, the principal factors deter- 
mining whether or not the material at the crater lip will be ejected are the 
pressure gradient tangential to the wave front and the shear stresses in the 
wave front. 

It is obvious that the choice of a cube as the beam or ring element on which 
the forces act is arbitrary. At the very high pressure gradients established 
near the impact point these forces will act on a very thin ring since the material 
nearer the impact point will have been ejected (Figure 29a). As the intensity of 
the shock decreases, the thickness of the ring on which these forces act will 
increase (Figure 29b). But the increase in thickness is limited because the 
pressure decreases sharply toward the open crater (Figure 29c). This increase 
in thickness is consistent with the observation that small fragments are ejected 
from the crater early in the process at high velocities and large fragments at 
low velocities later in the process. Reasonable dimensions for the radial 

distance over which the pressure difference acts to cause shear are also in- 
dicated by the thickness of the crater lip after crater formation is complete 
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(Figure 29d). 


3. The Dynamic Shear Modulus and Intrinsic Strength of Structural 
Materials, An examination of the elastic property data in Tables 4, 5, and 7 
indicates that high values of dynamic shear modulus are observed for low den- 
sity metals such as beryllium, refractory metals such as tungsten, and for 
ceramics such as corundum (Al^O^) and silicon carbide which contain elements 
from the center of periods 2 and 3 of the periodic table. If one calculates the 
intrinsic shear strength as a fraction of the shear modulus as suggested by 
Frenkel and Mackenzie, the low density metals and ceramics will have the 
highest strength to specific gravity ratios. The intrinsic shear strength ( G/Zrf ) 
to specific gravity ratio of beryllium is 1.85 x 10^ psi and for silicon carbide 
it is 1 . 23 x 10^ psi. A comparable value for 34T stainless steel is 
0. 22 x 10^ psi. 

There remains some uncertainty with respect to the choice of shear 
modulus values. The values chosen here are for polycry stall ine samples. 
However, almost all single crystals show a substantial degree of elastic aniso- 
tropy. The effect of anisotropy on crater formation has been illustrated by 

(55) 

Eichelberger and Gehring. If the impacting particles are much larger than 

the individual grains of the target material so that the forces resulting from 
impact are spread over many grains, the correct choice would seem to be the 
properties of the polycrystalline body. In intermediate cases, in which failure 
can occur due to forces spread over only a few grains, the choice is uncertain. 
Some materials are very weak in shear in particular crystal directions. For 

/ c Z \ 

example, the shear moduli for beta brass (CuZn) given by Lazarus' are 

C 44 = 8.24 x 10 11 and 1/2 (C n - C 12 ) = 0.975 xlO 11 dynes /cm 2 . 

cm 
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This large anisotropy is characteristic of body centered cubic structures com- 

57 

posed of metals whose ions have closed shells (Zener ). In the case of large 
grain size and substantial anisotropy, it may be important to use the smallest 
shear modulus to calculate the intrinsic strength. 

4. Variation of the Dynamic Shear Modulus with Temperature and 
Pressure . The effect of temperature and pressure on the shear modulus, may 

be important in calculation of the intrinsic strength. The temperature is of 
direct importance in the case of space radiators operated at temperatures high 
enough to affect the mechanical properties of the materials used. The shear 
modulus of materials usually decreases with increasing temperature but in a 
few cases, for example vitreous silica, an increase is observed. The tempera 
ture variation of the shear modulus is influenced by the change in shear modului 
with volume at zero temperature, the thermal energy and the thermal expansior 
The effect of target temperature on crater volume has been illustrated by 

/ C C \ / CO \ 

Eichelberger and Gehring' ' based mainly on data given by Allison et al. ' 
The results indicate a marked dependence of crater volume on the properties 
of target materials as these properties are affected by temperature. The meta! 
used were lead, aluminum, copper, zinc, and cadmium. These targets were 
impacted at 5.01 Km/sec. The crater volumes increased substantially with 
increasing target temperature. 

The variation of the shear modulus of pure aluminum as a function of 
temperature has been studied by K§. 1 The shear modulus decreases linearly 

with temperature for aluminum single crystals and there is a pronounced de- 
crease in shear modulus of polycrystalline aluminum in the temperature range 
from 200-400°C. Internal friction measurements at low frequencies indicate 
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that this variation in polycrystalline aluminum is caused by grain boundary 
relaxation. 

The internal friction measurements indicate a long relaxation time for the 
grain boundary movement. Therefore, the linear variation in shear modulus 
shown for the single crystal sample is more appropriate than the values for 
pure polycry stalline aluminum for use with samples at the high strain rates 
characteristic of meteorite impact. 

Similar data for beryllium has not been found so far. For equiaxed grains 
of uniform size, Zener^ ^ has shown that the ratio of the relaxed to the unre- 
laxed modulus is given by 

__ 2j7±5j0 

in which j> is the Poisson's ratio. The Poisson's ratio of beryllium is very 
low and shows marked anisotropy, varying from 0.10+ 0. 02 J."c" axis to 
0. 035 + 0. 01 || "c" axis. Therefore, the ratio of the relaxed to unrelaxed shear 
modulus (Gj^/G^j) will be between 0.42 and 0.45. The variation of the dynamic 
elastic modulus (Young's modulus) with temperature is given by White and 
Burke. ^ ^ It decreases from about 45 x 10^ psi at room temperature to about 

39.5 x 10 psi at 1200°F. If the temperature variation of the shear modulus 
is similar to the temperature variation of Young's modulus, the combined effect 
of the ordinary temperature variation and grain boundary relaxation could lead 
to a relaxed shear modulus at 1200°F as small as 37% of the room temperature 
value. As mentioned earlier, the unrelaxed values are probably applicable at 
the high strain rates characteristic of meteoroid impact. The characteristics 
of the grain boundaries can be altered by alloying if this is considered necessary. 
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The pressure dependence of the shear modulus may also be important. 

The important pressure for the crater wall criterion is that existing in the 

material near the crater wall at the time the crater formation ceases. These 

pressures are relatively low compared with the pressures existing at earlier 

times. Since the pressures are lowest near the crater lip and the elastic 

constants usually increase with the increasing pressure, the weakest material 

subject to shear forces is the material at the crater lip. The pressure range 
4 5 

of interest is 10 to 10 bars. 

The pressure dependence of the elastic constants has been measured by 

several investigator s . Measurements have been made for KCi, 

NaCZ, CuZn, Cu, Ag, Au, Al, oc SiO^ and Mg; in most cases to pressures of 

about 10 kilobars, the lower end of the range of interest. 

If we take impacts of aluminum on copper as an example to evaluate the 

importance of the pressure dependence of the shear modulus, we find that the 

4 

shear modulus of copper can be expected to increase about 2.8% in 1 0 bars 

5 

and about 28% in 1 0 bars. The pressure derivative of the shear modulus of 
other materials is larger than that of copper in many cases but in these cases 
the shear modulus and critical shear strength are usually smaller so that the 
pressure range is smaller and the percentage change in the shear modulus for 
this pressure range may remain about the same. These changes in shear modu- 
lus are somewhat smaller than the uncertainties in the choice of a strength 
criterion. Furthermore, the increases due to pressure may be somewhat offset 
by decreases due to temperature rise when the shockwave passes. Therefore, 
the use of the shear modulus at atmospheric pressure in the critical shear strength 
criterion is justified until further improvements are made in this shear strength 
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criterion. The most important variation in properties likely to occur in space 
radiators is the decrease in the elastic constants due to high operating tempera- 
tures. For aluminum used at 370°C 700°F) the decrease in shear modulus 

for a single crystal is about 23%. In the case of steel impacting on 25-0 
aluminum, the increase in temperature from room temperature to 370°C 
results in a 19% increase in crater radius. Again, in view of the crudeness 
of the crater-wall criterion it doesn't seem necessary to make a correction 
but, at higher temperatures, a correction would probably be necessary. 
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III. SUMMARY AND ILLUSTRATIVE CALCULATIONS 

This part presents a brief synopsis of the present status of this theory, 
and shows how to use the formula for crater size. No familiarity with other 
parts of the report is assumed. The formulas below are presented in such a 
way that the interested reader may apply them readily, without the necessity for 
going through the lengthy analysis above. 

1 . Review of the Analysis 

When a particle strikes a target surface at high speed, large amounts of 
energy and momentum are quickly deposited over a very small portion of the 
surface. This release drives a strong shock wave, usually hemispherical in 
shape, into the target, generating extremely large pressures. Because these 
pressures (normal stresses) are so large compared with the material’s resist- 
ance to shear deformation, one may consider that the target material behaves 
like an inviscid, compressible fluid. Thus the problem of determining the 
response of the target is essentially that of solving for the flow behind a strong, 
hemispherical shock wave. This flow pattern is described by the standard 
fluid-mechanical equations expressing the conservation of mass, momentum, 
and energy, together with the equation of state of the medium. The solution of 
these equations is an exceedingly formidable task. They are coupled, nonlinear 
partial differential equations in two independent spatial variables, as well as the 
time. Furthermore, as the shock wave advances into the target, distributing 
its energy over more and more mass, the strength of the shock decays, so that 
ultimately the fluid-mechanical approximation itself ceases to be valid. At this 
stage, a transition to a model which properly accounts for the plastic and elastic 
behavior of the medium is needed. 
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In order to make possible an analytic solution of this complicated prob- 
lem, a series of approximations to the physical picture has been made. They 
can be grouped into three main areas, and are described in the paragraphs that 
follow. 

The first simplification is to examine only that class of solutions in which 
the fluid motion is self- similar, i. e. , the distribution of physical quantities at 
successive instants is assumed to be the same when viewed on a scale defined 
by the shock radius at each instant. The mathematical consequence of this is 
to suppress time as an independent variable. Each physical quantity, instead 
of depending separately on the time and on the distance r from the impact 
point, is now a function of the combination r/R s (t) , where R s (t) is the 
shock -wave radius at any instant. This reduction in the number of independent 
variable constitutes a significant simplification in the differential equations that 
must be solved. At the same time, it imposes two restrictions. 

The first restriction, which assures that there will be no explicit time de- 
pendence in the differential equations, is that the shock radius must grow as a 
power of the time. Thus, out of all the possible solutions of the basic equations, 
we examine only that class for which R s Cn t* • 

The second condition imposed by the similarity assumption is that the 
boundary conditions must also be independent of the time, and this in turn re- 
stricts the type of state equation that can be used to those of the form 

e = 

where e is the internal energy per unit mass, -fV the pressure, and p is any 
function of the density p . In actual fact, the thermodynamic states of typical 
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target materials are well represented by the Mie-Griineisen equation of state, 
in the range of interest here. This equation may be written in the form 


e _ pKp) A ^ 

where the Gruneisen constant, P , depends weakly on , and where the 
function A depends only on density. The leading term of this equation can be 
accommodated in a self -similar solution, and will obviously be a good approx- 
imation at very high pressure, where the second term, A.(/0) > is negligible. 

The theory in its present state assumes that A(^>} can be neglected, and 
makes the further approximation of using the equation of state of a perfect gas 
of constant specific -heat ratio 7 

'f' 

e = \)f> 

This approximation amounts to the neglect of the small variation in the 
Gruneisen factor F(^o) > and permits all the known analyses for blast waves 
in perfect gases to be applied to the present problem. In every impact, some 
portion of the hydrodynamic flow will occur at pressures sufficiently low that 
the influence of the term A(f) is not negligible. In an attempt to account par- 
tially for the presence of the term A(^) • the value of that is used in any 
particular impact case is taken as the value that matches the actual state equa- 
tion at the impact point. For a perfect gas, the density ratio across a strong 
shock is always /) . Thus, if the actual density ratio across the 

shock that starts into the target from the impact point were to be 2, then T' 
would be chosen as 3, For a given projectile -target combination, the "effective 
•/" depends on the impact speed. 
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The second important area in which approximations have been made deals 
with the extent to which the conservation of energy and momentum are satisfied. 
Two types of solution have been studied. The first, which allows for spatial 
variations in two directions, permits conservation of energy and momentum, and 
explicitly describes the pattern of mass ejection from the target. The second 
type of solution considers only variations radially from the impact point, con- 
serves only the total energy, and does not describe the pattern of mass ejection. 
All the results that have been found for the former solution indicate that it is 
nearly identical to the vastly simpler constant-energy solution. The physical 
explanation for this is not presently known, although several speculations have 
been advanced. The simplest of these is based on the experimental observation 
that targets struck by hypervelocity projectiles often acquire momenta many 
times that of the projectile, implying that the material ejected from the target 
must also carry several times the projectile momentum. Thus it appears the 
momentum of the projectile itself plays a relatively minor role in the overall 
conservation process. 

These first two approximations -- that of similarity, and that dealing with 
neglect of momentum conservation -- are sufficient to permit an extremely 
simple solution for the time history of the shock radius 


R s (t) = 


f 25 
{dn'I l (r) 



/s 


where E is the kinetic energy of the projectile, ± the time after impact, 
the undisturbed target density, and a function representing the influence 

of the equation of state. This prediction shows reasonable agreement with the 
very limited amount of data that have been published so far. It disagrees chiefly 
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, which 


. 2/3 

in that the experiments tend to favor a time dependence more like Z 
suggests that the experiments lie in a range intermediate between the strong- 
shock limit ( R s ^ t^ 5 ) and the acoustic limit ( t ). In any event, the 

agreement between experiment and the present form of the theory is not bad, and 
may be expected to improve at higher impact speed. Thus, if used with caution, 
it will serve as a basis for crater-size predictions. 

This brings us to our third and final approximation, which is not asso- 
ciated with the fluid-mechanical equations as such, but rather with the later 
stages of crater formation, where the plastic and elastic response of the target 
becomes important. This stage is reached when the inviscid-fluid approxima- 
tion breaks down, i. e. , when the target’s resistance to shear deformation is no 
longer negligible compared to the pressure being generated by the shock. Thus, 
the blast-wave solution must be cut off at the instant when this pressure has 
become equal to the intrinsic shear strength, Gr/ZTf , & being the dynamic 

shear modulus. Ideally, one would use the conditions at the instant of cutoff as 
the starting values for a plastic- and elastic -wave solution, from which the 
ultimate crater radius could be found, as a function of the cutoff value of the 


shock radius, among other things. The theory in its present form does not ac- 
tually carry out such a solution. Instead, the radius of the crater is taken to 
be equal to the cutoff value of the shock radius. The formula that results 





is quite simple, and agrees well with experiment, especially in materials like 
lead, in which hypervelocity conditions are achieved at relatively low impact 
speeds. 
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In summary, then, the net effect of the set of approximations described 
above is to make possible a simple solution, giving the dependence of the crater 
radius on the pertinent impact parameters. In addition to its utility in making 
predictions, the analysis contributes to our fund of analytical understanding of 
the phenomenon as a whole, especially as regards the importance of momentum 
conservation, and the separate influences of the state equation and material 
strength. The theory is capable of considerable improvement, by accounting for 
nonsimilar effects in the fluid-dynamic phase, and by accomplishing a better 
description of the plastic and elastic stages of the process. The need for more 
definitive experiments on shock-wave time histories is clearly indicated, and, 
furthermore, it appears that some measurements of the dynamic strength 
properties of targets must be made if the results of experiments are to be inter- 
preted properly. 

2. Illustrative Calculations 

The formula recommended for crater size is Eq. (78) of Section I. 

Rc _ ( 1 V 

ci \/2(7+l)I'(/) F 

Here denotes the crater radius (the crater is assumed to be hemispherical) 

and d is the diameter of the projectile. (For nonspherical projectiles, d 
should be taken as the diameter of a sphere of equal mass. ) is the density 

of the projectile material, and V the impact speed. The symbol IP denotes 
the intrinsic shear strength of the target, Gy/ 27 ?' , where &■ is the dynamic 

shear modulus as measured, for example, by an ultrasonic technique. Table 4 
lists values of G-/<?7 f for a number of materials. Both Pp^ and ^ have 
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the units of pressure. A convenient set of units to work with is to take f> p in 

3 2 2 

gm/cm , V in cm/sec, and IP in gm/cm-sec , i. e. , in dynes /cm . Some 

useful conversion factors are 


Multiply : 
feet/ sec 
km/ sec 
psi 

megabars 


By 

30. 48 
10 5 

6. 895 x 10 4 


T o Get 
cm/sec 
cm/ sec 
dynes / cm^ 
dynes /cm^ 


The only quantity remaining to be determined is 7 , which characterizes the 
state equation of the target. The quantity 7 is actually the adiabatic index of 
a perfect gas whose equation of state matches that of the target in the vicinity of 
the impact point. The determination of 7 is thus essentially a solution for the 
initial conditions that exist between the shock wave that propagates into the tar- 
get and the shock wave that travels back into the projectile. The steps to follow 
are given below 

1. Choose a value for the pressure in the region between the shocks. 

This pressure is the same in the target and projectile. 

2. Find the particle speeds in both the target and projectile materials by 
applying the following equation to each one, respectively 



2 , 

Table I gives the values of C , S , and p o c for various materials. 

3. The impact speed, V , is now found by adding the two velocities 

found in the preceding step: V= U-,) + , ) • This is the speed 

e V 'PROJECTILE ''TARGET 

required to generate, at the impact point, the pressure chosen in Step 1. 
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4. The value of Y corresponding to the impact speed of Step 3 is now 
found from 

r= 2 & * sm j _ , 

wherein the values used are those appropriate to the target material. 

Repeating these steps at a series of values of then determines T'(V) , 
and this, in turn,allows one to calculate » using Fig. 8 for the function 

(r+l) I, (r) . 
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CONCLUDING REMARKS 


The net result of the set of approximations described above is to 
make possible an extremely simple solution, giving the dependence of 
crater radius on the diameter, density, and speed of the projectile, and 
on the dynamic strength and state equation of the target. Both fluid- 
mechanical and solid-state considerations play a role in this solution, and 
their individual effects are combined by the simple expedient of terminating 
the blast-wave solution at the point where the shock pressure has decayed to 
the intrinsic strength of the target. The fact that the solution for the fluid- 
mechanical phase is simple enough to permit such an incorporation of the 
target strength is a direct consequence of the fact that a constant-energy 
approximation can be used. The formula that results displays reasonable 
agreement with experiment, and may be expected to improve at higher impact 
speeds, where the strong- shock approximations hold over a longer portion 
of the process . 

Improvements in both the fluid and solid aspects of the theory describ ed 
above can be achieved, at little or no expense in simplicity of the results. 

In regard to the blast-wave portion, it would appear that the most fruitful 
area to explore is the development of non- similar solutions. The most 
definitive test to which this part of the solution can be put is the prediction 
of shock location as a function of time. The experimental evidence presently 
available in this area is quite limited, and apparently lies in a transitional 
regime where the shock speed has begun to decay from the blast-wave behavior 
*-> t 2|/5 ) toward the acoustic limit (/? s • To properly account for this 
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behavior, the hydrodynamic phase of the solution must be extended, so as 
to incorporate the nonsimilar effects that enter by way of the low-pressure 
term in the state equation. In spite of this deficiency, however, even the 
presently available form of the solution gives a surprisingly accurate pre- 
diction of the shock-wave location, and it may be expected to improve at 
higher impact speeds. 

On the solid-state side of the question, the most obvious area in need 
of improvement is the rule for determining the crater size. There is no 
theoretical obstacle to making the transition from blast-wave theory to a 
plastic- and elastic-wave description at late time. Such an analysis would 
provide a general relation between the final crater dimensions and the shock 
radius at the instant of transition. Presently, these are assumed equal. 

Finally, there is a need for further study of the importance of 
momentum conservation, since the neglect of this consideration is crucial 
in achieving the simplicity of the present solution. 

Certain aspects of these theoretical findings have application to the 
design and interpretation of experiments. Quite apart from its obvious use 
as a practical method of estimating penetration, the implications of this 
theory are that many experimental conditions must be carefully examined 
if results of experiments are to be properly interpreted. In particular, it is 
important to know the dynamic strength of the target material. Furthermore, 
the influence of the equation of state demonstrated in these studies produces, 
at high speed, a variation of crater radius with impact velocity suggestive of 
"momentum scaling". The fact is that the total energy is the dominant para- 
meter throughout the range of impact speeds, while reduction of the crater 
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size below the value appropriate to pure "energy scaling" can be attributed 
to the effect of the equation of state. Thus, extensive experiments in sub- 
stances whose state equations are not well known would appear to be incapable 
of determining the relative contributions from these two sources. 
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APPENDIX A 


TRANSFORMATION OF THE BLAST -WAVE EQUATIONS 


(A-l ) 


This appendix outlines the steps by which the centerline equations 
may be transformed to the , 2 coordinate system, where: 

. _ N $ , = Nr - 3(0 f n 

*~~f- ’ ' 1 " 1 se ' 

For the constant-energy solution, of course, Y , and hence 2 , are zero. 

A convenient starting point is to use the centerline equations in the form 

jl\ 

(A-2) 


+ Vd <t> +• '/'(% * 2 - 0 

-irN. = 0 

-2-^- 

The quantities <£> and -f can be eliminated from these, using 

d<f> = dr i) 


d f 


(A- 3) 
(A-4) 


4L .Ax. + zA± + 

f * y * 1 t 

The resulting set of equations contains only ^ ^ and 2- • 

(fL -N) dL&n> f* +■ aLsC + (3d + 2 ) cL&v fj ~ 0 


(A- 5) 


(A-6) 


(A- 7) 


-t+fr-N) + d +zr-r — a- +Z — — = 0 


(A -8) 
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(A-9) 


— 2 . y + (fiC — Ni) 




Equation (A-7) has been used in deriving the form of (A-9) shown here. If 

aLin, & 

we now use the continuity equation (A-7) to eliminate — from (A-8), 

and the energy equation (A-9) to eliminate , a little algebra 


leads to 


-%-(3r*+r* fr-n)fr-o 

^^7 ('■ t-Nf-'f 

which may be regarded as a means of finding f| , once # , y 

I dy- 

are known. Next, since the quantity — can be written as 

I dy _ cCJtwy d# 
y ~d7JUbfi «Ld cCdm/f^ 

equation (A-9) can be rearranged to 

<d£ _ -& y-Q * 

r- / ~+ (g- a/) 

Equating this value of dd^dXmf! to that given by (A-10) then gives 

I dy __ \fc-N f- y\ - Q- 1) a, -(3y- 1 ) y- / 

JT d * (d-N) 

In these coordinates, the equation for dcd/d$ may be written as 


(A-10 ) 


and H 


(A-l 1 ) 


(A-l 2) 


where 


_ z ~"j[ g2 ~ 


%e ~ 


_N* 5>V 

if 1 f 99* 




(A-l 3) 
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APPENDIX B 


SIMULATION OF METEOROID IMPACT BY ENERGY RELEASE 


A major conclusion reached in Part I, above, is that crater formation 
is controlled chiefly by the energy of the impacting particle, its momentum 
playing only a secondary role. Thus we may expect to simulate hypervelocity 
impact by any experiment in which a strong shock wave is driven into a 
target by the deposition of energy in any form. 

It is of central importance, in considering any simulation of this type, 
to be certain that the mode of energy deposition does in fact drive a strong 
shock wave into the target. We shall return to this question below, but for 
the moment we assume that this condition has been achieved, and present the 
formulas that follow as a consequence. 

The relation between crater radius and the amount of energy absorbed 
by the target is Eq. (77): 


R 


c 


I 


jf(r+l)X { (r) 



(B-l ) 


To specify Y in the particle-impact case, we matched the actual conditions 
behind the shock at the impact point, and found that Y depended only on the 
impact speed V .We may follow the same procedure for the case of energy 
deposition by some other means, except that we must now identify the 
experimental parameter that fixes the shock strength, and hence Y . The 
quantity that does this is the power being absorbed by the target, per unit 
area in the plane of the shock. To see why this is so, consider a plane shock 
wave of unit area advancing at speed U s into a medium of undisturbed 


RM- 1 655-M-4 


112 



density p . 


P, » «, 

* o 




(°o ’ Ao ’ e o = o 

^o=0 


In unit time, this shock processes an amount of mass given by p o m s per 
unit area, and raises it to the energy (per unit mass) 


V* \ Pi t 


(B-2) 


Thus the rate of energy acquisition by the material behind the shock, per unit 
time and area, is 


power /area = f> 0 U . , 


H-a-)- 


T I*' u 


(B-3) 


5 

The strength of any shock wave may therefore be characterized by the amount 

of power per unit area which it delivers to the medium through which it 

travels. In particular, such a specification of shock strength serves to define 

y . The Hugoniot curve for iron is interpreted in this light in Fig. 22, where 

it is seen that weak shockwaves, characterized by Y - 10 (i.e. , = 11/9) 

9 2 

impart about 10 watts/cm , while extremely strong shocks, for which Y - 2, 

13 2 

transfer to the medium some 10 watts/cm . These orders of magnitude 
are typical of metals. In fact, for any material whose Hugoniot is given by 


U s - C + S u, 


(B-4) 


a little algebra will show that 

3 I 


power/area = p C 


2 lr-(2s-/> 



l + 


2s 


Y- (2s -I), 


(B-5) 
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This function is shown in Fig. 23. In view of the fact that 1 . 0 megabar-km/ sec 

is equal to 10^ watts /cm^, it can be seen that shock waves in the range of 

9 13 2 

interest here will generally have "power ratings" from 10 to 10 watts/cm . 

It is interesting, in this connection, to note that the experiments reported by 

9 11 , 

Altshuler et al achieved shock waves of strength equivalent to 4 x 1 0 watts /cm 

If we consider the deposition of a total energy E over an area A of the 
target surface in a time Y , it is clear that f will be a function only of E/At* . 
Thus the craters formed, according to Eq. (B-l), will depend on the two para- 
meters E and Ar • Fig. 24 shows the craters predicted by Eq. (B-l) for 
deposition of energy in iron targets, at various values of E and AY . Over- 
laid as dotted lines are analogous predictions for the craters formed by impact 
of iron projectiles of various sizes, striking at various speeds. As noted above, 
there is a one-to-one correspondence between impact speed and the power 
density at the impact point, and both of these are shown on the abscissa scale 

of this figure. The point to be noted is that any experimental technique capable 

11 2 

of driving shockwaves of strength greater than 10 watts /cm can simulate 

impact conditions which are presently beyond the capability of conventional 

projection techniques. One energy source that immediately suggests itself for 

such an application is the laser. By focussing the beam from such a device, 

13 2 

power densities of 1 0 watts/cm , delivered in less than a microsecond, can 
be achieved^ with existing equipment. The fact that the maximum output of 
these devices is currently being improved at such a rapid rate indicates that 
simulation by a laser beam is a promising experimental technique. 

The calculations presented in Fig. 24 to illustrate the simulation 
possibilities start from the hypothesis that the energy absorption takes place 


RM-1 655-M-4 


114 



by means of a blast wave mechanism. Particularly in the case of electro- 
magnetic energy deposition, this assumption needs careful scrutiny. There 
would appear to be little doubt that is is the correct mechanism when the 
rate of energy input is sufficiently high. Certainly, the different phenomena 
brought into play, at one extreme, by striking a match and, at the other, 
by detonating a nuclear bomb, provide evidence that the mechanism of energy 
absorption changes, at some point, from one of linear heat conduction to the 
nonlinear shock-wave mechanism. Exactly where such a transition will 
occur in solid media is not presently known, although it is presumably 
amenable to theoretical analysis. The conclusions reached above are based 

on the assumption that a shock wave will be the correct mechanism when- 

11 2 

ever the incident power density exceeds 10 watts /cm . 
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APPENDIX C 


SOLUTION NEAR THE SONIC POINT 


At the "sonic" point, two quantities must be zero: 

(**-N 


(C-1) 


(C-2) 


The asterisk is used here and in the equations that follow to denote conditions 
at the sonic point. Because these quantities are zero, the values of 
and are indeterminate, according to Eqs. (A-12) and (A-10). To 

resolve the indeterminacy, put 


7 Z 7 


(/+€); ^ (l +&) ) 2 - = 


and let 


K d fs 


f) ® A > (: 


d± \ = a 

cU.) ~ ^2 


so that 


*-■*>%* s ’ & 


(C-3) 


(C-4) 


(C-5) 


. y/ 

While all three of 4 * , y and appear in these equations, it should 

be noted that only two of them are independent, since (C-1) and (C-2) must 
hold simultaneously. 




By using Eq. (C-3), the expression for 

■=-(/-!)+ 




may be written as 








- 

)**- 2 ( N + i)**+n] 


(C-6) 
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A second relation between and may be found from 

V* 


o _ fdj \ _ {oi ) 

2 "^/ ' (d*/dt~ n )* 





1 ^ (#*-N) i 


V. 

t*-Nf (3*Jb 2 ) \ 

*% n;j 


(C-7) 


Thus, (C-6) may be written as 


(C-8) 


which gives a linear relation between ^ and sk z . Using this relation in 
(C-6) then gives a quadratic for Jk t : 


where 


A/^, + + C - 0 

(C-9) 

a=f 5 + F„ 

(C-l 0) 

8 - F 3 - -F i F s +F , _F i F» 

(C-ll) 

C = -F,F, + P,F t -F 2 F s 

(C-12) 

F, = -(r-0(ic*- N) 

(C-l 3) 

f 2 = (V'-m) [2 -(/-/) ?“-(3 /) * *] 

(C-14) 

F 3 = 2fc*-N) 

(C-l 5) 

F 5 = 3-*:*Vz* - ^-(l— N) 

(C-l 6) 

F t - 3**-2(N+l)«*+N 

(C-l 7) 
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F. 




**('-¥ -**)-* tie 

z-(r-i)z*-(3r-i)** 


2 -(r-l)i* -3(r-i)** 


Having and ^2 , one can now find 

(**-Nf(3+JkJ + A,(3 **+*''-i g-N]) - Fi 
"U t)J 2 



(**-n) 31 . 
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INTERNAL ENER6Y, (KILOMETERS/ SEC) 
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Figure 4 HUGONIOT CURVES FOR IRON AND LEAD, AFTER ALTSHULER ET AL 
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Figure 5 EFFECTIVE 2* FOR LEAD TARGETS STRUCK BY IRON 
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Figure 6 ATTEMPTED SPHERICALLY-SYMMETRIC SOLUTION, N =0.30 
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Figure 7c DISTRIBUTION OF DENSITY, PARTICLE VELOCITY, AND PRESSURE 
CONSTANT-ENERGY, SPHERICALLY-SYMMETRIC BLAST WAVES 

7 * = Y 
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PROJECTILE DENSITY 

SHOCK PRESSURE AT WHICH THE BUST-WAVE 
SOLUTION IS CUT OFF 
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Figure 10 CONSTANT-ENERGY, EFFECTIVE- 7 * PREDICTION OF CRATER SIZE IN SEMI-INFINITE TARGETS 
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P = 0.2x10® psi (0/30) 
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EFFECTIVE T CHOSEN SO AS TO NATCH CONDITIONS AT THE IMPACT POINT 4 ..; 
EXPERIMENTAL DATA FROM HERRMANN AND JONES, REF. 69 AND 70 \ j j i | 


IMPACT VELOCITY V , km/ttc 

Figure 13 CONSTANT-ENERGY, EFFECTIVE- 3* SOLUTION FOR 
CRATERS FORMED BY ALUMINUM STRIKING COPPER 
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Figure 14a THE FUNCTION ^(V) FOR FUSED QUARTZ 
STRIKING VARIOUS TARGETS 
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Figure IHb THE FUNCTION 3*(V) FOR FUSED QUARTZ 
STRIKING VARIOUS TARGETS 
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Figure me THE FUNCTION ^(V) FOR IRON STRIKING VARIOUS TARGETS 



Figure md THE FUNCTION ^(V) FOR IRON STRIKING VARIOUS TARGETS 
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V, km/sec 


Figure I4f THE FUNCTION T'(V) FOR ALUMINUM 
STRIKING VARIOUS TARGETS 
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SOLUTION FOR N = 0.37S BASED ON THE APPROXIMATION 



Figure 17 CONSTANT- ENERGY, SPHERI CALLY-SYMMETRI C SOLUTION (N = 0.4) COMPARED 
WITH "CENTERLINE" SOLUTION (N = 0.375) 
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Figure 19 DEPENDENCE OF N ON 7FOR ONE-DIMENSIONAL BLAST WAVES 
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Figure 23 POWER-DENSITY RATING FOR MA 
IS DESCRIBED BY U=C+S> 



CT SPEED SHOWN 
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Figure 24 BLAST-WAVE PREDICTION OF CRATERS FORMED IN IRON BY DEPOSITION 
OF ENERGY E OVER AREA A IN TIME X 
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Figure 25 HUG0NI0TS FOR NORMAL AND POROUS IRON 
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megabars 












PRESSURE 



AGRAM OF POROUS 


(b) LATER TIME 



(c) p vs DISTANCE FROM THE IMPACT POINT (AT LATE TIME) 



Figure 29 EJECTION OF MATERIAL FROM A CRATER 
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Table I 
HUGONIOT DATA 
U $ = C + Su f 

I ^71 ^ I 1 ja I 0/2 77 

MATERIAL gm/cm 3 km/ sec S MEGABARS MEGABARS | SOURCE 

SILVER 10.49 3.22 1.64 1.094 

GOLD 19.24 3.06 1.61 1.800 

CADMIUM 8.64 2.41 1.72 .501 

COBALT 8.82 4.65 1.51 1.909 

CHROMIUM 7.10 5.18 1.54 1.902 

COPPER 8.90 3.97 1.48 1.404 0.069 

MOLYBDENUM 10.20 5.17 1.20 2.73 

NICKEL 8.86 4.67 1.41 1.930 

LEAD 11.34 2.07 1.52 .484 0.0069 

TIN 7.28 2.67 1.43 .518 

THORIUM 11.68 2.08 1.38 .505 

TITANIUM 4.51 4.79 1.07 1.033 

THALLIUM 11.84 1.82 1.57 .393 

VANADIUM 6.1 5.11 1.21 1.592 

TUNGSTEN 19.17 4.00 1.27 3.07 0.21 

ZINC 7.14 3.04 1.58 .661 

BISMUTH 9.80 2.0 1.36 .392 

IRON 7.87 4.0 1.59 1.259 0.14 

ANTIMONY 6.62 2.0 1.67 .265 

BERYLLIUM 1.82 7.98 1.09 1.157 0.21 

MAGNESIUM 1.74 4.49 1.27 .351 

INDIUM 7.31 2.37 1.61 .411 0.0062 

NIOBIUM 8.57 4.45 1.21 1.695 0.062 

PALLADIUM 12.0 3.79 1.92 1.726 

PLATINUM 21.45 3.67 1.40 2.89 

RHODIUM 12.44 4.68 1.64 2.72 

TANTALUM 16.6 3.37 1.15 1.890 

ZIRCONIUM 6.5 3.77 0.93 .924 

ALUMINUM 2.7 5.85 1 . 1 1 .924 0.041 

LUCITE 1.18 2.59 1.51 .0792 0.0021 

FUSED 2.20 1.30 1.56 .0372 0.048 

QUARTZ 
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Table 2 

N(^) FOR ONE-DIMENSIONAL BLAST WAVES 


7 

r + 1 

N 

SOURCE 

7- l 

1 

OO 

.5 

REF. 28 

i.i 

21 

.56888 

If 

1.4 

6 

.6 

H 

5/3 

4 

.61073 

m 

2.8 

2. II II 

.626704 

it 

CO 

1 

.6416 

n 

I. OOI 

2001 

.5106 

PRESEHT REPORT 

I.I 

21 

.5683 

II 

1.2 

1 1 

.5843 

II 

1.3 

7.6667 

.5935 

N 

1.4 

6 

.6000 

If 

1.8 

3.5 

.6143 

II 

2.0 

3 

.6182 

n 

2.5 

2.3333 

.6244 

fi 

3.0 

2 

.6279 

n 

4.0 

1.6667 

.6321 

ii 

5.0 

1 .5 

.6341 

it 

7.0 

1.3333 

.6365 

N 

10.0 

1 .2222 

.6385 

II 

100.0 

1 .0202 

.6412 

If 

1000.0 

1 .0020 

.6416 

If 
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Table 3 

STAGES OF ENERGY TRANSFORMATION DURING IMPACT 


STAGE 

ENERGY TRANSFORMATION 

DURATION OR TERMINATION 

»• ELASTIC DEFORMATION 

K.E. TO STORED MECHANICAL ENERGY 

ENDS WHEN INTERFACE PRESSURE 
EXCEEDS DYNAMIC YIELD 
STRENGTH 

2. PLASTIC DEFORMATION 

K.E. TO ABSORBED MECHANICAL 
ENERGY 

ENDS WHEN INTERNAL ENERGY 
EXCEEDS LATENT HEAT OF 
SOLIDIFICATION 

3. SOLID STATE PHASE 
CHANGE 

K.E. TO ABSORBED LATTICE ENERGY 

OCCURS IN SOME CRYSTALLINE 
SOLIDS AT ELEVATED TEMP. 

& PRESSURE 

1. HIGH PRESSURE 

ABSORBED MECHANICAL ENERGY TO 

CONTINUES UNTIL THERMAL 

COMPACTION 

THERMAL ENERGY 

ENERGY EXCEEDS REMAINING 
KINETIC ENERGY 

5. SOLID/LIQUID AND 

ABSORBED MECHANICAL ENERGY TO 

OCCURS WHEN LOCAL TEMPER- 

LIQUID/ GASEOUS 
PHASE CHANGES 

THERMAL ENERGY 

ATURE EXCEEDS HIGH PRESSURE 
MELTING OR BOILING POINT 

6. GASEOUS EXPANSION 

ABSORBED THERMAL ENERGY TO 

OCCURS WHEN LOCAL RATE OF 

(EXPLOSION) 

MECHANICAL ENERGY 

ENERGY ABSORPTION EXCEEDS 
RATE OF ENERGY RELEASE BY 
MECHANICAL OR THERMAL RAD- 
IATION 

7. BLAST WAVE 

ABSORBED THERMAL ENERGY TO 
ABSORBED AND STORED MECHANICAL 
ENERGY PLUS MECHANICAL ENERGY 
RADIATED AWAY FROM TARGET 

EFFECTIVELY ENDS WHEN RATE 
OF ENERGY RELEASE FALLS 
BELOW THE STRONG SHOCK WAVE 
VELOCITY 

8. GASEOUS/ LI QUID AND 

ABSORBED THERMAL ENERGY TO 

OCCURS AS THE OUTER TEMPER- 

LIQUID/ SOLID 
PHASE CHANGE 

MECHANICAL ENERGY 

ATURES OF THE ENLARGED 
DAMAGE AREA FALL BELOW THE 
BOILING AND FREEZING POINTS 

9. PLASTIC DEFORMATION 

DYNAMIC MECHANICAL ENERGY TO 
ABSORBED MECHANICAL ENERGY 

OCCURS WHEN BLAST WAVE 
PRESSURE EXCEEDS DYNAMIC 
YIELD STRENGTH. ENDS WHEN 
STRAIN RATE BECOMES LESS 
THAN WEAK SHOCK WAVE 
VELOCITY 

10. SHOCK WAVES AND 
STRESS WAVES 

RADIATION OF MECHANICAL ENERGY 

CONTINUES UNTIL ALL DYNAMIC 
MECHANICAL ENERGY IS 
ABSORBED 

II. FRAGMENTATION 

STORED MECHANICAL ENERGY TO 
ABSORBED MECHANICAL ENERGY 

OCCURS AFTER RAPIDLY DEVEL- 
OPED HIGH COMPRESSIVE 
STRESSES 
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Table 4 

PHYSICAL PROPERTIES OF SELECTED SOLIDS (AT ROOM TEMPERATURE) 


MATERIAL 

DENSITY 

f 

gm / cm^ 

MELTING OR 
SOFTENING 
POINT 

°C 

YOUNG'S 

MODULUS 

Y 

I0 6 psi 

BULK 

D 1 LATATIONAL 
WAVE VELOCITY 
C l 

km / sec 

DYNAMIC 

SHEAR 

MODULUS 

G 

1 0 6 psi 

SHEAR 

WAVE 

VELOCITY 

C S 

km / sec 

INTRINSIC 

SHEAR 

STRENGTH 

1 TO i- 
30 2X 

I0 6 psi 

PLASTICS 

POLYETHYLENE 

0.90 

IOO°(S.P.) 

0. 1 1 

2.0 

0.055 

0.5 

0.002-. 009 

POLYSTYRENE 

1.06 

86°(S.P.) 

0.76 

2.4 

0.19 

1 .1 

0.006*. 03 

UICITE* 

1.18 

85°( S.P. ) 

0.58 

2.7 

0.21 

l.l 

0.007-.03 

SOFT METALS 
INDIUM 

7.28 

156° 

1 .57 

2.4 

0.54 

0.7 

0.02 -.09 

LEAD 

11.3 

327° 

2.3 

2.0 

0.78 | 

0.7 

0.03 -.1 

S0LDER(6OSn-9OPb) 

8.52 

188° 

4.5 

3.1 

1.58 

l.l 

0.05 -.3 

METALS 

ALUMINUM 

2.71 

660° 

10.5 

6.4 

3.6 

3.0 

0.1 - .6 

OORAl ( I7ST) 

2.67 

590° 

10.4 

6.3 

3.9 

3.1 

0.1 - .6 

COPPER 

8.9 

1083° 

18 

5.0 

6.7 

2.3 

0.2 -1 

STAINLESS STEEL (39T) 

7.9 


28.4 

5.8 

1 1.0 

3.1 

0.4 -2 

IRON 

7.7 

1535° 

29.9 

5.9 

11.9 

3.2 

0.4 -2 

TUNGSTEN 

19.3 

3910° 

52.5 

5.4 

19.4 

2.6 

0.6 -3 

BERYLLIUM 

1 .84 

1278° 

44.6 

12.9 

21 .3 

8.9 

0.7 -3 

GLASS 

PYREX 

2.23 

820° (S.P.) 

9.0 

5.6 

3.6 

3.3 

0.1 - .6 

FUSED SILICA (Si0 2 ) 

2.20 

I670°(S.P.) 

10.5 

5.9 

4.5 

3.7 

0.2 - .7 

CERAMICS 

HEMATITE (Fe 2 0j) 

5.20 

1565° 

34 

6.8 

15 

4.6 

0.5 -2 

MAGNETITE (Fe^) 

5.18 

1538° 

33 

7.3 

13 

4.2 

0.4 -2 

CORUNDUM (AIjOj) 

3.98 

2020° 

59 

II. 0 

23.5 

6.4 

0.8 -4 

SILICON CARBIOE 

3.1 

2600° 

58 

12.0 

24 

7.3 

0.8 -4 

POROUS SOLIDS 
IRON (P = 99.3*) 

0.05 

1535° 

_ 

- 

0.01 

- 

<0.001 

IRON (P = 93.5*) 

0.5 

1635° 

-< 0.3 

< 2 

0.1 

-C 1 

<0.01 

IRON OXIDE ( P=99 . 1 * ) 

0.05 

~ 1560° 

- 

- 

-c 0.01 

- 

<0.001 

(P=90.9*) 

0.5 

✓*. 1 550° 

< 0.3 

< 2 

-< 0.16 

-s 1 

<0.02 


* LUCITE IS A TRADE NAME FOR POIYMETHYL-METHACRYLATE 
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Table 5 

VARIATION OF ELASTIC CONSTANTS WITH PRESSURE 


MATERIAL 

ELASTIC 

MODULUS 

M 

PRESSURE COEFFICIENT 
OF ELASTIC MODULUS 

PRESSURE 

RANGE 

STUDIED 

REFERENCE 

0C-QUARTZ 

®I4 

-9.0 xO ' 6 cm 2 /k g 

1000 ATMOS. 

67 

(Si0 2 ) 

®44 

4.9 x 0 " 6 cm 2 /kg 




C 66 

- 6 . 8 x 10 “® cm^lkq 



FUSEO SILICA 

C ll 

-1.3x10'® cm^/dyne 

3.5x10® dynes/c * 2 

66 

(Si0 2 ) 

C 44 

- 1 . 0 x 10 “® cw^/dyne 



ALUMINUM 

K =i ( c ,| + 2C |2 ) 

0.5UI05 4^ ( b ; r ) 


46 (p.324) 


C 44 

0.75x10® ^ ( b -|>) 



COPPER 

« =i (C| , + 2C| 2 ) 

0.41x10® typ- ( b f r ) 


46 (p.324) 


€44 

0-31*10® ^-(b-ir) 
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Table 7 

COMPRESSIVE STRENGTH OF MATERIALS UNDER HYDROSTATIC PRESSURE 

(AFTER BRIDGEMAN) 52 
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Table 8 

INTRINSIC STRENGTH BASED UPON SHOCK WAVE DATA 
(BASED UPON McQUEEN AND MARSH 54 ) 


MATERIAL 

C CRITICAL PRESSURE 
( p H = (Tt) AT dP/dV = 0 

SHEAR 
STRESS 
& s max 

DYNAMIC 
SHEAR 
MODULUS Q 

s max 
G 

Ag 

-174 

kbar 

2.52xJ0 6 ps i 

I.26xl0 6 psi 



Au 

292 

H 

¥.23 

ii n 

2.12 

II n 



Be 

273 

ii 

3.96 

n ii 

1.98 

n ii 

8.9 

1/4.5 

Cd 

75 

ii 

1.09 

n it 

0.55 

ii ti 



Co 

371 

ii 

5.43 

ii n 

2.72 

ii ii 



Cr 

330 

n 

¥.79 

it n 

2. ¥0 

n n 



Cu 

233 

n 

3.38 

ii n 

1.69 

ii ii 

6.7 

i/4 

Hg 

32 

ii 

0.46 

n it 

0.23 

ii ii 



Mg 

70 

n 

1.01 

n ii 

0.51 

n ii 



Mo 

548 

H 

7.95 

ii n 

3.98 

n ti 



Nt 

331 

n 

4.80 

ii it 

2.40 

ii n 

11.6 

1/4.8 

Pd 

77 

ii 

I..I2 

it ii 

0.56 

ti n 

0.78 

1/1.4 

Sn 

86 

it 

1.25 

n ii 

0.63 

ii n 



Ti 

237 

n 

3.44 

n n 

1.72 

n ii 



T 

68 

» j 

0.99 

ii n 

0.50 

ii it 



V 

329 

n 

4.77 

n ii 

2.39 

ii ii 



W 

606 

it 

8.80 

it ii 

4.40 

it it 

19. ¥ 

1/4.4 

Zn 

107 

ii 

1.55 

ii ii 

0.78 

n it 
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